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ABSTRACT 

Heat- transfer  measurements  were  made  for  a  single, 
water-heated  tube  in  a  pool  of  R-114  to  simulate  operating 
conditions  of  a  water  chiller.  Data  were  obtained  for  a 
smooth  copper  tube,  and  for  two  commercially  available 
tubes:  a  spirally  roped  copper-nickel  tube  with  a  porous 
coating;  and  a  copper  tube  with  a  structured  outer  surface 
and  a  multiple-start  helical  ridged  inner  surface. 
Measurements  were  made  for  refrigerant-oil  mixtures  at  oil 
concentrations  from  0  to  6  mass  percent  with  a  boiling  pool 
temperature  of  13.  8  C.  Results  for  the  two  enhanced  tubes 
with  and  without  oil  are  compared  to  the  smooth  tube  data. 
Enhancement  factors  for  the  overall  heat- transfer  coeffi- 
cient were  4.  0  and  3.  6  for  the  structured  surface  and 
porous-coated  tubes,  respectively,  in  pure  refrigerant  and 
at  a  water  Velocity  of  2  m/s.  For  these  same  conditions  the 
enhancement  factors  for  the  outside  heat- transfer  coeffi- 
cient were  14.6  and  6.4  for  the  porous-coated  and  structured 
surface  tubes ,  respectively. 
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I.  INTRODUCTION 

Recent  developments  in  boiling  surfaces  have  shown 
considerable  enhancements  in  heat  transfer  performance.  The 
worldwide  literature  on  enhanced  heat  transfer  contains  over 
3000  published  technical  papers  and  reports  [Ref.  1:  p.  81] . 
This  increased  interest  in  heat  transfer  augmentation  is  the 
result  of  incentives  for  energy  and  material  savings.  One 
effective  way  to  improve  the  heat  transfer  is  by  using 
passive  augmentation. 

Passive  augmentation  uses  fine-scale  alteration  of 
surfaces,  both  external  and  internal  [Ref.  1:  p.  82].  These 
surfaces  may  consist  of  either  an  applied  porous  coating  or 
fins  deformed  in  various  ways  to  provide  a  large  number  of 
reentrant  cavities  on  external  surfaces.  Methods  of 
internal  enhancement  include  insert  devices,  forged  fins,  or 
deformation  of  the  surface  (i.e. ,  spirally  roped  or  corru- 
gated tube).  The  use  of  such  surfaces  can  lead  to  consider- 
able reduction  in  the  size  and  weight  of  heat- transfer 
equipment.  The  reductions  may  lead  to  smaller  capital  costs 
or  operating  costs  or  both. 

Of  particular  interest,  in  naval  applications,  is  the 
reduction  in  the  size  of  water  chillers  in  refrigeration 
systems.  While  a  number  of  investigations  are  currently  in 
progress,  the  theoretical  treatment  of  the  boiling  perform- 
ance of  various  tubes  is  almost  impossible  owing  to  the  very 
complicated  mechanisms  involved.  The  degree  of  difficulty 
increases  with  the  presence  of  oil  in  the  refrigerant 
liquid.  In  general,  refrigeration  systems  with  oil- 
lubricated,  hermetically  sealed  compressors  contain  small 
mass  percents  of  oil  in  the  refrigerant  liquid.  Therefore, 
reliable  data  covering  a  wide  variety  of  operating  condi- 
tions  for   various   refrigerants   and    different   boiling 
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surfaces  is   essential  to  further   develop  more   compact  and 
efficient  evaporators. 

Most  of  the  experiments  on  enhanced  boiling  surfaces 
reported  in  the  literature  have  used  electrically  heated 
tubes.  The  difficulties  associated  with  these  tubes,  espe- 
cially during  the  instrumentation  stage,  raise  questions  as 
to  the  reliability  of  the  data  as  described  by  Wanniarchchi 
et  al  [Ref .  2:  p.  14] .  Despite  the  precautions  taken  to 
minimize  contact  resistance,  it  may  still  be  present  between 
the  thermocouple  locations  and  the  outer  boiling  surface. 
Another  area  of  doubt,  described  by  Wanniarchchi  et  al. 
[Ref.  2:  p.  14]  ,  is  the  nonuniform  heat  fluxes  generated  by 
commercially  available  electrical  heaters.  Additionally,  in 
order  to  use  an  electrical  heated  tube,  the  inside  surface 
of  the  tube  must  be  smooth,  or  any  internal  enhancement  must 
be  bored  out.  This  is  disadvantageous  as  designers  would  be 
interested  in  knowing  the  overall  heat-transfer  coefficient. 
Most  industrial  experiments  use  2- to  3-meter- long  tubes 
with  warm  water  flowing  inside  them.  While  the  information 
generated  from  these  experiments  closely  approximate  actual 
evaporators ,  compared  to  electrical  heated  tubes ,  the  large 
water  temperature  drop  from  inlet  to  exit  (up  to  5  K)  may 
make  it  difficult  to  study  the  details  of  the  boiling 
process. 

The  type  of  refrigerant  selected  will  also  influence  the 
design  of  evaporators.  R-114,  a  moderate-pressure  refrige- 
rant, is  receiving  more  attention,  in  particular  for  naval 
applications.  Advantages  to  using  R-114  include:  (a)  it 
belongs  to  the  refrigerant  group  with  the  lowest  toxicity, 
(b)  it  is  very  stable  with  temperature,  and  (c)  it  has  a 
fairly  large  value  for  the  energy  transfer  per  unit  volume 
of  vapor  (i.e.,  Ev  =  hfg  /vg  )  [Ref.  3:  pp.  281-291].  The 
advantages  of  the  first  two  items  are  clear.  The  third  item 
is   advantageous  since   a  higher   value  for   Ev  means   lower 
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pressure  drop  through  the  refrigerant  piping  for  a  specified 
heat  duty.  As  the  normal  boiling  point  of  the  type  of 
refrigerant  increases,  the  Ev  value  increases.  R-ll,  R-113, 
R-114,  and  R-22  have  Ev  values  in  increasing  order.  R-114 
may  be  preferable  to  R-ll  and  R-113  due  to  the  higher  Ev 
value  and  preferable  to  R-22  due  to  the  ability  to  use 
lighter  components  in  the  refrigeration  loop. 

Based  on  the  above  discussion,  the  major  objective  of 
the  present  investigation  was  to  study  the  boiling  perform- 
ance of  two  commercially  available  tubes,  a  porous-coated 
tube  and  a  mechanically  structured  surface,  in  comparison 
with  a  smooth  tube  for  the  following  conditions:  (a)  tubes 
are  water-heated,  (b)  boiling  fluid  is  R-114  with  0,  1,  2, 
and  6  percent  by  mass  oil,  and  (c)  boiling  temperature  is 
13.8  °C. 


13 


II.     BACKGROUND 

A.  NUCLEATE  BOILING 

When  the  temperature  of  a  boiling  surface  exceeds  the 
saturation  temperature  of  a  fluid  by  a  few  degrees,  nucleate 
boiling  occurs.  The  difference  between  the  boiling  surface 
and  saturation  temperature  is  the  amount  of  wall  superheat. 
Dougherty  [Ref.  4:  p.  175]  defined  pool-boiling  as  vaporiza- 
tion occurring  under  the  following  conditions:  (1)  liquid 
depth  >>  the  bubble  diameter,  (2)  there  is  negligible  effect 
on  heat  transfer  due  to  low  or  no  externally  imposed 
velocity  of  the  fluid,  and  (3)  the  bubbles  move  away  from 
the  boiling  surface  due  to  a  force  field.  With  the  above 
conditions,  the  heat  flow  can  be  expressed  as: 

h  -  -/fa  (2-D 

A  description  of  the  boiling  mechanism  is  given  by 
Chongrungreong  and  Sauer  [Ref.  5:  p.  701].  A  thin  layer  of 
superheated  liquid  is  formed  adjacent  to  the  boiling 
surface.  Bubbles  nucleate  in  this  thin  layer  and  grow  from 
preferred  spots  on  the  boiling  surface.  It  is  assumed  that 
the  primary  resistance  to  heat  transfer  is  within  the  thin 
liquid  layer.  The  height  of  the  liquid  in  the  boiling 
container  is  not  a  primary  variable,  but  the  liquid  proper- 
ties should  be  controlled.  However,  large  flooded  evapora- 
tors may  experience  a  significant  "submergence  effect"  from 
the  liquid  head. 

B.  EXTERNAL  SURFACE  ENHANCEMENT 

For  smooth  heating  surfaces,  bubbles  nucleate  at  various 
scratches  and  cavities  on  the  surface.  Fujii  stated  that  the 
number  of  active  sites  increases   as  the  heat  flux  increases 
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[Ref .  6:  p.  48]  .  Webb  [Ref .  7:  p.  46]  stated  that  the 
ability  to  enhance  the  nucleate  boiling  coefficient  by 
applying  some  type  of  roughness  has  been  known  for  over  50 
years. 

Presently,  there  are  two  major  types  of  commercially 
developed  roughnesses:  (1)  porous  coating  and  (2)  various 
fins  and  surfaces  with  reentrant  cavities  [Ref.  8:  p.  24] . 
The  first  type  of  roughness  is  a  sintered  particle  coating, 
usually  copper  or  aluminum.  Figure  2. 1  shows  a  schematic  and 
a  photomicrograph  of  a  porous  coating.  The  small  reentrant 
cavities  are  interconnected  by  substrate  tunnels.  The 
surface  reduces  the  required  superheat  for  vapor  generation 
by  entrapping  a  high  density  of  relatively  large  vapor 
nuclei  in  the  cavities  contained  within  the  porous  coating 
[Ref.  8:  p.  24].  As  described  by  Webb  [Ref.  7:  p.  58], 
researchers  found  that  a  critical  pore  size  not  a  particle 
size  governs  the  number  of  nucleation  sites.  Large  pores  are 
required  for  liquids  with  high  surface  tension  and  high 
thermal  conductivity,  while  small  pores  work  best  for 
liquids  with  low  surface  tension  and  low  thermal  conduc- 
tivity (i.e. ,  refrigerants).  Czikk  and  O'Neill  [Ref.  9:  p. 
53]  developed  correlations  to  include  the  effects  of 
different  pore  sizes  of  porous-coatings.  They  concluded  that 
there  were  two  resistances  to  heat  transfer  for  these  coat- 
ings: a  nucleation  superheat  related  to  bubble  diameter;  and 
a  conduction  superheat  controlled  by  the  liquid  film  sepa- 
rating the  vapor  bubbles  from  the  metal  matrix.  Carnavos 
[Ref.  10:  pp.  106-108]  found  that  a  porous  coating  resulted 
in  a  700-800  percent  better  heat- transfer  coefficient  than  a 
smooth  surface  in  a  pool  of  R-ll.  Czikk  [Ref.  11:  p.  98] 
also  found  increased  performance  of  the  porous-coated  tubes 
used  in  his  20-ton  water-chiller  experiment. 

Structured  fins   or  surfaces  form  reentrant   cavities  of 
various  geometries.    Figure  2. 2   is  shows   an  example   of  a 
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Figure  2.  1. 


Schematic  and  Scanning  Electron  Micrograph 
(50Qx)  of  High  Flux  Surlace  (irom  Ret.  3). 
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Reentrant  Angle 


Figure  2.2.  Schematic  of  Manufactured  Reentrant  Cavity 
C from  Ref.  3 ). 
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manufactured  reentrant  cavity.  These  reentrant  cavities  act 
as  stable  nucleation  sites,  which  enhances  the  heat 
transfer.  The  nucleation  site,  to  remain  active,  is  depen- 
dent on  the  mouth  diameter  falling  within  a  critical  range 
and  shape  with  a  maximum  reentrant  angle.  This  range  is  a 
function  of  the  fluid  properties.  The  cost  of  these 
surfaces,  presently,  are  not  significantly  higher  than  the 
cost  for  smooth  tubes  and  they  dramatically  improve  boiling 
performance  compared  to  the  smooth  tubes. 

The  mechanisms  by  which  the  reentrant  cavities  operate 
are  complicated.  Previous  research  concluded  that  the  combi- 
nation of  bubble  evaporation,  thermal  boundary- layer  strip- 
ping and  bubble  agitation  controlled  the  heat  transfer  from 
smooth  surfaces  [Ref.  12:  p.  192].  In  contrast,  different 
experiments  conducted  by  Arshad  [Ref.  12:  p.  192]  and  Arai 
[Ref.  13:  p.  37]  show  that  thin  film  evaporation  in  the 
reentrant  cavity  controlled  the  heat- transfer  mechanism. 
Bubbles  were  formed  by  vapor  exiting  cavities  as  the  liquid- 
vapor  interface  of  the  thin  film  contacted  the  cavity 
surface.  Surface  tension  holds  most  of  the  liquid  on  the 
cavity  walls.  Ayub  [Ref.  14:  p.  64]  showed  this  similar 
"thermosiphon  mechanism"  with  his  experiments  using  enhanced 
surfaces. 

Nakayama  [Ref.  15:  p.  37]  gives  a  fairly  detailed 
description  of  a  physical  model  undergoing  this  bubble 
growth  mechanism.  Figure  2. 3  shows  the  three  major  phases  of 
bubble  growth.  Phase  I  consists  of  a  pressure  buildup  in 
the  tunnel  by  evaporation  of  liquid  held  in  corners  and 
continues  until  the  meniscus  reaches  a  hemispherical  shape 
of  radius  r0  =  d0/2  (where  d0  is  the  mouth  diameter  of  the 
reentrant  cavity).  In  Phase  II,  the  meniscus  grows  faster  at 
some  pores  than  at  others.  The  Phase  I  pressure  buildup  is 
reduced  as  the  vapor  enters  the  growing  bubbles.  The 
meniscus  at  inactive   pores  does  not  grow  due   to  this  vapor 
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Figure  2. 3.  Physical  Model  of  Bubble  Mechanism 
(  from  Kef.  15). 
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pressure  reduction.  Initially,  the  bubble  expands  under  the 
high  internal  pressure,  while  lacer  expansion  is  controlled 
by  the  receding  liquid  inertia  around  the  bubble.  Phase  III 
is  termed  the  liquid  intake  phase  with  liquid  flowing  into 
the  tunnel  through  inactive  pores.  This  flow  is  over  a  short 
interval  of  pressure  depression  which  occurs  as  the  pressure 
of  the  bubble  and  tunnel  is  lower  than  the  pool  pressure. 
The  bubble  leaves  the  pore  and  new  meniscus  formation  closes 
off  the  pore,  ending  this  phase  and  returning  the  cycle  to 
Phase  I. 


C.   EPFECTS  OF  OIL 

The  introduction  of  oil  into  the  pure  refrigerant,  in 
general,  causes  a  decrease  in  the  boiling  coefficient  as  the 
oil  concentration  increases.  The  effects  of  oil  on  the 
mixture  are  theoretically  complicated  as  most  empirical 
equations  relating  physical  properties  are  derived  form 
experimental  data  and  not  from  theoretical  equations. 
Various  empirical  equations  are  cited  by  Chongrungreong  and 
Sauer  [Ref.  5:  pp.  703-705].  Yet,  Jensen  and  Jackman 
[Ref.  16:  p.  184]  showed  that  the  correlation  developed  by 
Chongrungreong  and  Sauer  appeared  to  overpredict  as  oil 
concentrations  increased  due  to  poor  prediction  of  the 
mixture  viscosity. 

Predictive  equations  were  developed  by  Jensen  and 
Jackman  [Ref.  16:  pp.  136-137]  for  density,  viscosity, 
surface  tension  and  specific  heat  for  pure  R-113,  pure  oil, 
and  R-113-oil  mixtures.  These  equations  are: 
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o  =  a    +  (o  -o  )CC  *S  (2.  4) 

m    r     o   r   c 


c^  =  Ci-C  )c   ♦  c  (2.5) 

rm      c   ~rl   bol 

These  equations  are  based  on  a  physical  model  shown  in 
Figure  2.4  developed  by  Jensen  and  Jackman  [Ref.  16:  pp. 
187-188].  The  model  shows  a  bubble  growing  on  a  heated 
surface  as  the  refrigerant  evaporates  from  the  superheated 
liquid  phase  into  the  bubble  interior.  Due  to  the  less  vola- 
tile nature  of  the  oil  in  the  mixture,  it  does  not  evaporate 
into  the  bubble,  but,  remains  at  the  liquid-vapor  interface. 
The  decrease  in  the  rate  of  bubble  growth  is  due  to  the 
decrease  of  oil  diffusion  into  the  liquid  mixture  and  the 
decrease  of  refrigerant  through  the  oil  layer  into  the 
bubble.  The  decrease  of  the  bubble  growth  rate  decreases 
the  heat- transfer  rate. 

D.   INTERNAL  ENHANCEMENT 

Internal  enhancement  increases  the  heat- transfer  area, 
creates  additional  turbulence,  and  produces  secondary  flows, 
all  of  which  contribute  to  the  increase  in  the  heat 
transfer.  The  laminar  sublayer  is  assumed  to  be  the  prin- 
cipal resistance  to  heat  transfer  [Ref.  17:  p.  6].  The 
developed  secondary  flow  thins  the  sublayer  and  moves  it 
into  the  stream  where  turbulence  prevails  and  without  an 
increase  in  shear;  likewise,  heat  transfer  is  increased 
without  an  increase  in  friction  [Ref.  18:  p.  30] .  Most 
studies  show  that  the  friction  factor  does  increase  with 
internal  enhancement  but  the  amount  is  strongly  dependent  on 
the  internal  enhancement  geometry.  The  increase  of  both  the 
friction  and  heat  transfer  is  dependent  on  fin  or  rib  pitch, 
groove  depths,  and  helix  angles  [Refs.  17,19:  pp.  6,  19-21]. 
Additionally,   these   studies   have   shown   that  there   are 
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critical  geometrical  dimensions  that  give  the  maximum  heat 
transfer.  For  example,  short  fins  (<<  tube  diameter)  with 
rifling  will  raise  heat-transfer  coefficients  with  no  prob- 
lems of  flow  stagnation.  Yet,  medium  fins  may  develop  flow 
stagnation  and  decrease  heat- transfer  coefficients  [Ref.  19: 
pp.  22-23] .  As  such,  the  selection  of  internal  enhancement 
will  generally  improve  the  performance  of  boiling  surfaces. 
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III.  EXPERIMENTAL  APPARATUS 

A.   OVERALL  APPARATUS 

A  schematic  of  the  experimental  apparatus  is  given  in 
Figure  3.  1.  The  photographs  in  Figure  3. 2  show  two 
different  views  of  the  apparatus.  Complete  details  on  the 
design  and  construction  of  the  refrigerant  and  oil  compo- 
nents are  given  by  Karasbun  [Ref .  20:  pp.  24-32] .  A  discus- 
sion of  the  modifications  for  the  water  heating  mode  is 
provided  in  this  paper.  The  major  components  of  the  appa- 
ratus are:  two  Pyrex  glass  tees,  an  R-114  liquid  reservoir, 
a  water-ethylene  glycol  mixture  sump,  an  R-12  refrigeration 
system,  a  vacuum  pump,  a  water  supply  tank,  three  centri- 
fugal pumps,  a  flow  meter,  three  heaters,  a  graduated  oil 
cylinder,  and  an  oil  reservoir.  The  R-114  boiling  and 
condensation  occurred  in  the  lower  and  upper  glass  tees, 
respectively.  The  R-114  vapor  was  condensed  by  fhe  water- 
ethylene  glycol  mixture  pumped  through  a  copper  condenser 
coil  located  in  the  upper  glass  tee.  A  1/2-Ton  R-12  refrig- 
eration system  maintained  the  water-ethylene  glycol  mixture 
between  -18  and  -14   C. 

The  major  difference  between  the  experimental  apparatus 
used  by  Karasabun  and  Reilly  [Refs.  20,8:  pp.  24-32,  30-35] , 
and  the  one  used  during  the  present  investigation  concerned 
the  use  of  water  heating  instead  of  electric  heating  as  a 
means  to  provide  heat  to  the  boiling  tube.  Filtered  tap 
water  supplied  from  a  storage  tank  (15),  was  pumped  by  two 
1/2-hp  motors  (13  and  14)  connected  in  series.  A  three-way 
valve  (V-16)  was  provided  on  the  downstream  side  of  the 
second  pump  to  obtain  two  flow  paths.  The  normal  flow  path 
for  data  runs  was  through:  the  metering  valve  (V-17),  the 
flowmeter  (6),  the  inlet  mixing  chamber  (16),  the  boiling 
tube,  the  outlet  mixing  chamber  (16),   the  chambers  provided 
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with  heaters  (7  and  8),  and  returning  to  the  storage  tank 
through  valve  (V-20).  The  by-pass  flow  path  was  the  same  up 
to  the  three-way  valve  then  the  path  was  through  heater 
chamber  (8)  and  back  to  the  storage  tank  through  valve 
(V-20).  Valves  V-17  and  V-21  were  closed  in  order  to 
isolate  the  boiling  tube.  This  by-pass  option  was  essential 
for  preheating  the  water  to  higher  temperatures,  to  achieve 
higher  heat  flux  values  in  the  boiling  tube  compared  to 
values  obtained  from  room- temperature  water.  A  500-W  and  two 
1000-W  heaters  were  used  to  maintain  a  steady  water  inlet 
temperature  and  to  preheat  the  water.  Flexible  pressure  hose 
was  used  for  water  piping  throughout  the  apparatus.  The 
mixing  chambers  and  sections  of  hose  connecting  the  mixing 
chambers  to  the  boiling  tube  were  double  insulated  to  ensure 
accurate  temperature  measurements. 

Copper-constantan  thermocouples  were  used  to  measure 
thermal  emf ' s  for  R-114  liquid  and  vapor,  water  inlet  and 
outlet,  and  water-ethylene  glycol  mixture.  A  series- 
connected  thermopile  consisting  of  10  junctions  on  either 
end  was  used  to  measure  the  temperature  drop  of  the  water 
from  the  inlet  to  the  outlet  of  the  boiling  tube.  This  ther- 
mopile was  calibrated  against  two  quartz  thermometer  probes 
and  agreement  was  found  to  be  better  than  ±  0. 02  K. 

B.   BOILING  TUBE  CONSTRUCTION 

Figure  3.3  shows  schematics  of  the  smooth  tube,  and 
two  enhanced  tubes  used.  The  tubes  tested  were: 

1.  a  smooth  copper  tube; 

2.  a  90:10  copper-nickel  corrugated  tube  (commercially 
referred  to  as  Korodense  tube)  with  an  external, 
sintered  porous  coating  (i.e.  High  Flux); 

3.  an   internally   and  externally   enhanced   tube   (alloy 
C12200)    produced   by  Wolverine   Division   of   U. 0. P. 
( commerically  referred  to  as  Turbo-B); 

4.  a  porous-coated  (High  Flux)  Korodense  tube  with  the 
porous  coating  machined  off; 

5.  a  Turbo-B  tube  with  the  external  enhancement  machined 
off. 
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These  tubes  will  be  referred  to  as  High  Flux,  Turbo-B,  modi- 
fied High  Flux,  and  modified  Turbo-B  tubes,  respectively, 
with  the  exception  of  tubes  2  and  4.  Tubes  2  and  4  will  be 
referred  to  as  Korodense  tubes  with  respect  to  discussions 
of  only  the  internal  enhancement.  All  of  the  tubes  had  an 
active  (i.e.  ,  heated)  section  of  304.8  mm  in  length.  The 
smooth,  High  Flux,  and  modified  High  Flux  tubes  had  an  outer 
diameter  of  15.  9  mm  and  an  inside  diameter  of  12. 7  mm.  In 
the  case  of  the  High  Flux  tube,  these  diameters  represent 
the  values  in  the  smooth  portion  of  the  tube  (i.e. ,  before 
performing  the  corrugation  process).  The  Turbo-B  tube  had 
an  external  structured  surface  and  a  multiple-start  helical 
ridged  inner  surface.  This  tube  had  an  outer  diameter  of 
16.  9  mm  at  the  base  of  the  external  structured  surface  and  a 
minimum  diameter  of  14. 5  mm  at  the  internal  ridge  tips.  The 
modified  Turbo-B  tube  had  identical  diameters  as  the  Turbo-B 
tube.  The  63.5-mm-long  inactive  sections  on  either  end  was 
insulated  by  Teflon  sleeves  of  1. 6  mm  wall  thickness  located 
inside  the  boiling  tubes,  as  shown  in  Figure  3.3.  The 
requirement  for  using  these  two  modified  tubes  is  explained 
further  in  Chapter  IV  (EXPERIMENTAL  PROCEDURES). 

C.   COMPUTER- CONTROLLED  DATA  ACQUISITION  AND  REDUCTION 

Hewlett-Packard  equipment  was  used  for  data  acquisi- 
tion, control,  and  analysis,  as  well  as  for  storing  of 
data.  An  HP-3497A  acquisition  unit  was  used  to  read  the 
thermocouple  and  thermopile  outputs.  Table  1  gives  the 
channel  designations  used. 

Information  entered  by  keyboard  to  an  HP-9826A  computer 
unit  prompted  and  controlled  the  data  acquisition  unit. 
Data  were  analyzed  and  stored  with  the  same  computer  unit.  A 
step-by-step  description  of  the  data-reduction  procedure  is 
given  in  Appendix  A,  along  with  a  printout  of  the  program 
(DRP4)  used. 
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TABLE  1 
CHANNEL  DESIGNATIONS 

Channel  Function 

0-1  R-114  liquid  thermal  emf's  T(l) 

and  T(2) 

2  R-114  vapor  thermal  emf  T(3) 

3  water-ethylene  glycol  mixture 
thermal  emf  T(  4 ) 

4  water  inlet  thermal  emf  T(5) 

5  water  outlet  thermal  emf  T(b) 
20  thermal  emf  difference  across 

water  inlet  and  outlet  T(7) 
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IV.  EXPERIMENTAL  PROCEDURES 

A.  FLOWMETER  CALIBRATION 

This  section  gives  a  brief  description  of  the  flowmeter 
calibration;  a  more  detailed  explaination  is  given  in 
Appendix  B.  A  Fischer  Porter  flowmeter  was  used  to  indicate 
percent  water  flow  through  the  boiling  tube.  Prior  to 
conducting  any  boiling  runs,  it  was  necessary  to  develop  a 
correlation  relating  the  flow  percent  to  the  mass  flow  rate 
and  water  velocity.  Additionally,  any  temperature  effect  on 
viscosity  would  be  incorporated  into  this  correlation  as  a 
correction  factor.  Calibration  runs  were  conducted  over  a 
temperature  range  of  19  °C  to  38  C  and  varying  flowmeter 
settings.  The  correlation  (equation  (B.l))  was  determined 
from  these  results.  The  correction  factor  for  the  viscosity 
temperature  dependency  was  determined  to  be  1. 0  with  an 
accuracy  of  ±  0.02  as  explained  in  Appendix  B.  In  fact,  it 
was  not  possible  to  find  a  systematic  trend  for  this  correc- 
tion factor  with  the  water  inlet  temperature.  The  stated 
accuracy  mainly  consists  of  the  uncertainties  involved  with 
the  visual  reading  of  the  flowmeter  settings. 

B.  BOILING  TUBE  AND  APPARATUS  PREPARATION 

The  external  surface  of  the  boiling  tube  was  cleaned 
with  a  Nitol  (2%  nitric  acid  and  98%  ethyl  alcohol)  solution 
and  both  the  external  and  internal  surfaces  were  rinsed  with 
acetone  and  were  air  dried  prior  to  installation  in  the 
apparatus.  A  vacuum  test  at  an  absolute  pressure  of  about  50 
mm  Hg  was  conducted  on  the  refrigerant  side  of  the  apparatus 
after  the  boiling  tube  was  installed.  If  no  signs  of  leaks 
were  evident  after  two  hours,  the  system  was  pressurized  to 
a  gage  pressure  of  about  250  mm  Hg  with  R-114  vapor  and  the 
apparatus  was   checked  for  leaks   with  an   automatic  Halogen 
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Leak  Detector.  After  fixing  any  leaks  that  may  have  been 
present,  R-114  liquid  was  transferred  from  the  reservoir  to 
a  pre-determined  level  in  the  boiling  tee.  This  level  gave  a 
free  surface  40  mm  above  the  centerline  of  the  boiling  tube. 
At  the  set  level,  the  mass  of  the  R-114  liquid  was  computed 
to  be  2. 48  kg.  A  desired  saturation  temperature  of  4  C  and 
the  resultant  system  pressure  were  obtained  by  turning  on 
the  R-12  refrigeration  system  and  varying  the  flow  of  cold 
(approximately  -16  C)  water-ethylene  glycol  mixture  through 
the  condenser  coil  with  valve  V-9.  At  this  saturation 
temperature  and  pressure,  each  boiling  tube  developed 
nucleation  sites.  The  saturation  temperature  of  4  C  was 
maintained  for  at  least  30  minutes  prior  to  initiating 
nucleate  boiling  on  the  complete  active  length  of  the 
boiling  tube. 

The  initiation  of  nucleate  boiling  could  be  done  with 
cold  water  at  19  °C  or  with  warm  water  at  25  C  through 
valve  V-17  or  with  steam  from  a  steam  generation  system 
through  valve  V-22.  Early  experimental  data  showed  slight 
differences  (up  to  10  percent)  in  the  computed  outside  heat- 
transfer  coefficient  depending  on  the  highest  heat  flux  at 
which  the  boiling  was  initiated  on  each  boiling  tube.  As 
will  be  discussed  in  Chapter  V  (RESULTS  AND  DISCUSSION),  the 
steam  intiation  method  provided  the  highest  initial  wall 
superheat  and  ensured  the  best  possible  repeatability  of  the 
data.  This  intiation  method  was  used  for  all  of  the  runs 
unless  otherwise  specified. 

Steam  initiation  was  conducted  by  introducing  steam 
through  valve  V-22  and  discharging  steam  drain  valve  V-18 
for  about  one  minute.  After  this  period,  the  three-way  valve 
V-16  was  shifted  from  a  neutral  position  to  the  normal  flow 
position.  Valves  V-22  and  V-18  were  closed  and  valve  V-21 
was  opened  simultaneously  as  the  three-way  valve  was  shifted 
to  the  normal  position.   Valve  V-17  was  opened  wide  prior  to 
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shifting  of  the  three-way  valve.  As  these  valves  were  shift- 
ed, warm  water   at  25  °C  was  pumped  through  the  boiling  tube 
at  the  maximum  water  velocity,  which  in  turn  gave  the  maxi- 
mum heat  flux  for  the  given  water  inlet  temperature. 

C.   BOILING  DATA  RUNS 

Once  the  nucleate  boiling  initiation  was  completed,  the 
saturation  temperature  was  raised  to  13.8  ?C  by  adjusting 
valve  V-9.  This  valve  was  also  used  to  maintain  this  satura- 
tion temperature  within  ±0. 05  K  throughout  a  boiling  run. 
The  water  inlet  temperature  was  maintained  to  within  ±0.  08  K 
by  using  the  heaters  or  introducing  filtered  cold  water, 
respectively.  The  combination  of  heaters  and  cold  water  was 
used  to  maintain  the  desired  water  inlet  temperature  during 
each  run.  Two  types  of  boiling  runs  were  conducted  at  each 
oil  mass  percent. 

The  first  type  of  run  was  conducted  after  holding  the 
desired  saturation  temperature,  water  inlet  temperature  and 
initial  flowmeter  setting  at  steady-state  conditions  for 
about  10  minutes.  After  this  period,  the  flowmeter  setting 
was  decreased  with  two  data  sets  being  taken  at  each  flow 
setting.  The  two  data  sets  were  to  show  repeatability  at 
each  setting.  The  period  between  pairs  of  data  sets  was 
about  5  to  8  minutes,  with  steady  state  for  each  flowmeter 
setting  maintained  for  at  least  2  minutes.  The  second  type 
of  boiling  run  consisted  of  maintaining  the  previously 
mentioned  initial  steady-state  conditions  at  a  selected  high 
water  inlet  temperature.  The  run  was  commenced  by  intro- 
ducing cold  water  into  the  storage  tank  which  lowered  the 
water  inlet  temperature.  While  holding  the  water  velocity 
constant,  six  to  seven  data  points  were  taken,  on  a  one  time 
pass,  at  about  every  0. 3  K  inlet  temperature  decrease.  When 
the  inlet  temperature  reached  a  selected  low  temperature, 
the  steady-state   conditions  of   the  low   inlet  temperature, 
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water  velocity,  and  saturation  temperature  of  13.8  C  were 
maintained  for  about  10  minutes.  Then,  the  water  inlet 
temperature  was  increased  using  the  heaters,  while  holding 
the  water  velocity  constant.  Again,  data  points  were  taken 
on  a  one  time  pass,  in  about  0. 3  K  intervals,  until  the 
inlet  temperature  reached  the  previously  used  high  inlet 
temperature. 

Two  complete  boiling  runs  of  the  first  type  mentioned 
above  were  performed  at  0%  and  2%  oil  concentrations  to 
further  demonstrate  repeatability.  Table  2  shows  the  percent 
difference  in  wall  superheat  (i.e.  ,  AT)  and  heat  flux  for 
each  boiling  tube  at  these  two  oil  percents.  The  higher 
percentages  for  the  High  Flux  and  Turbo-B  tubes  for  the  wall 
superheats  are  simply  because  of  the  low  wall  superheat 
values  of  these  two  tubes  when  compared  to  the  values  for 
the  smooth  tube.  At  these  low  superheat  values,  a  small 
difference  in  the  wall  superheat  leads  to  a  larger 
percentage  difference  than  that  which  occurs  at  higher  wall 
superheats.  The  variation  of  the  water  inlet  temperature  was 
not  considered  a  factor  in  this  percent  difference  as  the 
variation  was  less  than  ±  0.02  K  per  minute  during  these 
constant-inlet- temperature ,  decreasing-water-velocity  runs. 

TABLE  2 
BOILING  RUN  REPEATABILITY 

Tube       Oil%      AT        q 


Smooth  0  ±  97  ±207, 

High  Flux  0  ±301  ±127 

Turbo-B  0  ±30$  ±  47 

Smooth  2  ±10%  ±257 

High  Flux  2  ±107  ±  27 

Turbo-B  2  ±30%  ±  2% 

Note:  initiations  performed  with  steam 


Additionally,  data  runs  were  conducted  to  determine  the 
temperature  increase  across  the  boiling  tube  due  to  the 
pressure   drop   between   the   locations   of   the   thermopile 
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probes.  The  data  collected  showed  that  this  difference  was 
less  than  the  ±0. 02  K  accuracy  of  the  thermopile  and  was 
considered  negligible. 

Following  the  boiling  runs  in  pure  R-114,  oil  was 
introduced  to  the  boiling  tee  through  valve  V-l  from  the 
graduated  cylinder  (5)  shown  in  Figure  3.1.  After  the 
required  volume  of  oil  for  a  given  oil  mass  percent  was 
added,  water  at  the  maximum  velocity  was  pumped  through  the 
boiling  tube  to  promote  vigorous  boiling.  This  boiling 
ensured  good  mixing  of  the  R-114  and  oil  mixture.  Either 
another  boiling  run  was  conducted  or  the  system  was  shut 
down  to  prepare  for  another  steam  initiation. 

D.   DATA  ACQUISITION  AND  REDUCTION 

The  automatic  data  acquisition  system  was  prompted  to 
record  the  required  thermal  emf ' s ,  by  keyboard  inputs  to  the 
computer.  The  data  were  immediately  processed  and  printed 
out  on  a  hard-copy  printer.  A  step-by-step  procedure  of  data 
processing  and  reprocessing  is  given  in  Appendix  A.  Also, 
a  printout  of  the  data  reduction  program,  DRP4 ,  is  included 
in  this  appendix.  The  heat  flux,  overall  heat-transfer  coef- 
ficient, outside  boiling  heat- transfer  coefficient,  and  wall 
superheat  values  were  computed  based  on  the  outside  area 
expressed  by  the  diameter  if  the  enhancement  of  the  outer 
surface  was  removed.  Further,  in  order  to  account  for  the 
heat  conduction  through  the  inactive  end  sections,  a  correc- 
tion was  made  for  heat  flux  using  an  iterative  computation 
procedure  based  on  natural  convection  as  discussed  by 
Karasbun  [Ref.  20:  pp.  54-56]. 

When  a  boiling  run  was  completed,  the  data  were  repro- 
cessed by  computer  to  obtain  the  inside  and  outside  heat- 
transfer  coefficients  for  further  calculations.  A 
Sieder-Tate- type  constant  (Cj)  for  the  inside  coefficient 
for  each  of  the  externally  smooth  tubes  was  calculated  by 
using  a  modified  Wilson  plot  I  Ref.  22j  .  For  the  smooth  tube, 
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modified  High  Flux,  and  modified  Turbo-B  tubes,  the  outside 
coefficient  was  calculated  simultaneously  using  the  cor- 
relation (equation:  (C.3)  and  (C.4))  developed  by  Rohsenow 
[Ref.  23:  p.  969]  ,  and  the  modified  Wilson  plot.  The 
Rohsenow  correlation  is  based  on  externally  smooth  tubes, 
thus  the  requirement  for  the  modified  High  Flux  and  Turbo-B 
tubes.  A  further  explaination  of  the  modified  Wilson  plot 
is  given  in  Appendix  C.  The  Cj  values  found  for  the  exter- 
nally smooth  tubes  were  used  for  the  corresponding  exter- 
nally enhanced  tubes.  Knowing  the  inside  coefficient,  the 
outside  coefficient  for  the  boiling  tubes  were  computed  by 
subtracting  the  inside  and  wall  resistances  from  the  meas- 
ured overall  thermal  resistance  (equation:  (C.7)). 
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V.  RESULTS  AND  DISCUSSION 

A.   INSIDE  HEAT-TRANSFER  COEFFICIENT 

Results  were  obtained  for  the  smooth,  High  Flux,  and 
Turbo-B  boiling  tubes  at  a  boiling  temperature  of  13. 8  C 
and  at  oil  concentrations  of  0,  1,  2,  and  6  mass  percent. 
The  heat  flux  versus  wall  superheat  for  these  three  tubes  is 
shown  in  Figure  5. 1.  All  three  tubes  are  in  the  nucleate 
boiling  range  with  decreasing  heat  flux.  It  can  be  seen  that 
the  High  Flux  and  Turbo-B  tubes  outperform  the  smooth  tube 
throughout  the  tested  heat  flux  range.  Additionally,  the 
Rohsenow  correlation  for  a  smooth  tube  at  a  boiling  tempera- 
ture of  13.8  °C  with  the  experimentally  determined  (from  the 
modified  Wilson  plot)  Csf  coefficient  of  0.0060  is  shown  in 
this  figure. 

The  Rohsenow  correlation,  as  used  in  the  modified  Wilson 
plot,  has  an  exponent  of  r  =  1/3.  Yet,  this  correlation  is 
very  sensitive  to  the  r  value.  A  change  in  this  exponent  of 
three  percent  yielded  a  change  in  the  Csf  only  of  0.016 
percent  but  gave  a  change  of  17  percent  in  the 
Sieder-Tate- type  constant  inferred  from  the  modified  Wilson 
plot.  For  the  purpose  of  data  analysis,  the  value  or  r  =  1/3 
was  used,  while  knowing  that  an  additional  uncertainty  was 
being  introduced  by  this  exponent. 

For  each  externally  smooth  tube,  the  Sieder-Tate- type 
constant  was  computed  using  the  modified  Wilson  plot  based 
on  5  to  10  different  runs  in  pure  R-114.  The  computed 
values  for  the  smooth,  Korodense,  and  Turbo-B  tubes  were 
0.036  ±0.002,  0.066  ±  0.005,  and  0.077  ±  0.003,  respec- 
tively. For  the  smooth  tube,  the  experimentally  found 
Sieder-Tate- type  constant  is  33  percent  greater  known  than 
the  well-known  value  of  0. 027  for  long  internally  smooth 
tubes   with  fully   developed   flow.     This  larger   constant 
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appears  to  be  the  result  of  the  entrance  effects  of  the 
shorter  experimental  tube  and  the  uncertainty  introduced  by 
the  exponent  r  of  the  Rohsenow  correlation  as  discussed 
above.  The  minimum  entrance  length  for  a  fully  developed 
pipe  flow  is  given,  as  a  rule  of  thumb,  by  Incropera  and 
DeWitt  [Ref.  24:  p.  406]  as 


—  =  60 

D 


(5.1) 


for  Reynolds  numbers  greater  than  or  equal  to  10000.  Also, 
as  the  pipe  roughness  increases  the  minimum  entrance  length 
decreases. 

Withers  has  developed  Stanton  number  correlations  for 
tubes  with  single-helix  and  multiple  helix  ridging 
[Refs.  25,26:  pp.  52-56,44].  These  correlations  are  listed 
below.   For  single-helix  ridging, 


St  = 


ft  Pr'0-5  (Z)°- 


7.22 


!.) Re  i 


f\0.5"  0.1 


(5.2) 


'  '   2..Sin[r  ^fl  (5-3) 

where  r  and  m  are  determined  by  tube  dimensions,  and 


Y  = 


2.5 


n  P  {1} 


*    3.75 


(5.4) 


For  multiple-helix   ridging, 


if       _    -0.5 
•  t    =-lT     Pr 


fl#    L     /f    N  0  .136 


(5.5) 


Z.46*n[r  W_)m] 
r  and  m  again  are  functions  of  tube  dimensions, 


(5.6) 
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Y  =  -^.Sln[2(fJ  ♦  3.75]j  (5.7) 

and 

1 
81  =  S.ssfsf  3  f£\°-i36  (5.3) 


A  comparison  of  Stanton  numbers  was  made  for  the  modified 
Korodense  and  modified  Turbo-B  tubes  using  the  following 
equations  for  these  two  tubes: 

h.  D. 
St  =   -i 1 (5.9) 


k^  Re  Pr 


where 


(5.10) 


The  results  for  these  two  modified  tubes  were  compared  to  a 
copper  Korodense  (i.e.,  single-helix  ridging)  and  a 
Turbo-Chill  (i.e.,  multiple-helix  ridging)  tubes,  respec- 
tively. The  latter  two  tubes  tested  by  Withers  were  similar 
in  radial  dimensions  but  were  1.  5  meter  long.  The  predicted 
St  numbers  using  equations  (5.1)  -(5.8)  for  the  Withers' 
Korodense  and  Turbo-Chill  tubes  were  0.00177  and  0.00273, 
respectively.  The  experimental  St  numbers  using  equations  ( 5.9 ) 
and  (5.10)  for  the  modified  Korodense  and  modified  Turbo-3 
tube3  were  0.00242  and  0.00236,  respectively.  The  majority 
of  the  difference  between  the  predicted  and  experimental  St 
numbers  may  be  attributable  to  the  entrance  effects  and  the 
uncertainity  introduced  by  the  Rohsenow  correlation. 

3.   LIGHT  OF?  EFFECTS 

As  referred  to  in  Chapter  IV,  Figure  5.2  shows  the  light 
off  effects  for  the  smooth,  High  Flux  (i.e.,  porous- coated 
Korodense),  and  Turbo-B  tubes.  These  effects  are  most  prob- 
ably due  to  the  requirement  of  greater  heat  flux  to  initiate 
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the  nucleation  sites  than  the  heat  flux  required  to  maintain 
the  sites  once  they  are  activated.  As  the  boiling  tube  sits 
in  the  R-114  liquid,  prior  to  initiation,  reentrant  cavities 
or  the  nucleation  sites  flood  with  liquid.  The  initial 
formation  of  the  vapor  bubbles  and  thus  the  initiation  of 
the  nucleation  site  requires  the  input  of  a  greater  heat 
flux  or  higher  wall  superheat.  Figure  5. 3  shows  the  boiling 
tube  initiation  with  steam  for  the  smooth,  High  Flux,  and 
Turbo-B  tubes  [Refs.  21,7:  p.  478,  62].  As  shown  in  Figure 
5.3,  the  High  Flux  tube  is  more  sensitive  to  light  off 
effects.  This  sensitivity  is  mostly  due  to  the  lower  wall 
superheats  at  which  nucleation  occurs.  A  change  of  .5  of  a 
degree  is  more  pronounced  at  these  lower  wall  superheats 
than  a  similar  change  at  the  wall  superheats  for  the  Turbo-B 
and  smooth  tubes. 

C.   SMOOTH  TUBE 

Figure  5.4  shows  the  performance  of  the  smooth  tube  at 
0,  1,  2,  and  6  mass  percent  oil  concentration  at  13.8  C. 
Reilly's  [Ref.  8:  p.  65]  data  for  the  electrically  heated 
smooth  tube  at  a  boiling  temperature  of  -2. 2  °C  in  pure 
R-114  liquid  are  also  included  for  comparison.  The  differ- 
ences between  the  experimental  results  shown  and  Reilly's 
data  are  probably  due  to  physical  characteristics  introduced 
in  the  electrically  heated  tube  when  soldering  the  copper 
sleeves  (inside  of  which  the  electric  heaters  were  fitted) 
on  the  interior  of  the  tube.  As  discussed  by  Reilly  [Ref.  8: 
pp.  57-61] ,  the  contact  resistance  at  the  interface  between 
the  sleeve  and  the  inner  surface  of  the  tube  was  minimized 
by  tinning.  However,  the  tinning  process  may  not  have  been 
100  percent  successful,  thereby  introducing  unacceptable 
uncertainties  into  the  experimental  measurements.  Further 
Reilly  reported  considerable  variations  (up  to  3  K  at  a  heat 
flux  of  98  kW/m2)  in  the  measured  wall  temperatures.  Reilly 
attributed  this  observation  to  an  axially  non-uniform  heat 
flux  generated  by  the  cartridge  heater. 
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As  shown  in  Figure  5.4,  the  wall  superheat  increased  for 
1  percent  oil  concentration  and  decreased  slightly  for  2  and 
6  percent  oil.  This  decrease  in  the  wall  superheat  for  the 
higher  oil  percentage  is  possibly  due  to  the  enhanced  bubble 
formation  experienced  by  the  tube  because  of  the  foaming 
action  of  the  R-114  and  oil  mixture  decreasing  resistance  to 
heat  transfer  as  explained  by  Chongrungreong  [Ref.  5:  p. 
701]  ,  Nobukatsu  [Ref.  27:  p.  60]  ,  and  Chaddock  [Ref.  21:  p. 
474]  . 

D.   HIGH  FLUX  TUBE 

The  performance  of  the  High  Flux  tube  in  terms  of  heat 
flux  versus  wall  superheat  is  given  in  Figure  5.  5. 
Additionally,  this  figure  shows  Reilly's  data  [Ref.  8:  p. 
70] ,  for  the  electrically  heated  High  Flux  tube  in  pure 
R-114  at  a  6.  7  C  boiling  temperature.  The  presence  of  oil 
increased  the  wall  superheat,  at  a  heat  flux  of  40  kW/m2 ,  by 
a  factor  of  1.  22  for  1  percent  oil  and  a  factor  of  2. 12  for 
6  percent  oil.  The  number  of  nucleation  sites  reduced 
considerably  and  the  foaming  action  increased  significantly 
as  the  oil  mass  percent  increased.  These  actions  are  shown 
in  Figure  5. 6.  With  6  percent  oil  concentration,  it  was 
visually  observed  that  the  nucleation  sites  were  active  only 
in  the  corrugations,  as  shown  in  Figure  5.6  (c).  This 
congregation  of  sites  is  due,  in  part,  to  the  sparser  porous 
coating  on  the  high  points  of  the  tube  and  the  decreased 
wall  thickness  in  the  corrugations.  As  a  result  of  the 
corrugation  process  during  tube  manufacture,  the  wall  thick- 
ness at  the  corrugation  is  smaller  than  the  rest  of  the 
tube,  which  implies  lower  wall  resistance.  During  the 
sintering  process,  the  coating  appears  to  have  concentrated 
in  these  corrugations  giving  a  thicker  coating  of  copper 
particles  and  increased  number  of  nucleation  sites  compared 
to  the  rest  of  the  tube  surface. 
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The  increase  of  the  wall  superheat,  when  oil  is  present, 
is  possibly  the  result  of  increased  surface  tension  and 
viscosity.  The  increasing  surface  tension  and  viscosity  lead 
to  the  formation  of  an  oil-rich  layer  next  to  the  boiling 
surface.  This  oil-rich  layer  has  an  insulating  and  inhib- 
iting of  bubble  growth  effect  within  the  cavities.  The 
higher  liquid-vapor  surface  tension  and  surface  tension  in 
the  oil-rich  layer  requires  more  energy  for  the  formation 
and  growth  of  the  bubble.  Additionally,  diffusion  processes 
affect  the  bubble  formation.  Oil  concentrations  near  the 
boiling  surface  are  high  compared  to  the  concentration  in 
the  bulk  liquid  and  diffusion  occurs  from  the  cavity  into 
the  bulk  liquid.  This  diffusion  continues  until  an  equilib- 
rium point  is  reached.  If  oil-rich  layers  are  developing  on 
the  boiling  surfaces,  the  diffusion  from  the  cavities  to  the 
bulk  liquid  slows  down  aiding  the  further  build-up  of  the 
oil-rich  layer.  The  vapor  bubble  forming  in  the  cavity  also 
has  an  oil  vapor  content  which  is  undergoing  its  own  diffu- 
sion process.  As  the  refrigerant  evaporates  through  the 
oil-rich  layer  into  the  interior  of  the  bubble,  it  must 
overcome  large  diffusion  resistances  set  up  by  the  oil  vapor 
within  the  bubble.  The  increased  surface  tension,  viscosity, 
diffusion  resistance  and  the  decreasing  bubble  formation 
rate  contribute  to  the  decrease  in  the  heat- transfer  coeffi- 
cient and  the  increase  in  the  wall  superheat  at  a  given  heat 
flux.   [Refs.  5,21,28,29,30:  p.  702,  477-478,59,372,82] 

Also,  boiling  runs  were  conducted  holding  the  water  vel- 
ocity constant  and  allowing  the  water  inlet  temperature  to 
decrease  and  then  to  increase.  The  decreasing  and  increasing 
heat  flux  with  the  respective  change  in  the  water  inlet 
temperature  versus  wall  superheat  is  shown  in  Figure  5. 7. 
The  motion  of  the  hysteresis  (i.e. ,  clockwise  or  counter- 
clockwise) depended  on  the  temperature  starting  point.  If 
the  temperature  was  decreased  and  then  increased,  the  motion 
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was  clockwise.  When  the  temperature  was  increased  and  then 
decreased,  the  motion  was  counterclockwise.  This  directional 
change  is  related  to  the  amount  of  heat  flux  required  to 
maintain  active  nucleation  sites.  In  the  first  case  of  tem- 
perature decreasing  then  increasing,  less  heat  flux  is  need- 
ed to  maintain  the  same  amount  of  sites,  but,  at  the  same 
time,  heat  flux  decrease  as  the  water  heat  capacaity  is 
decreasing.  In  the  second  case,  with  the  water  inlet  temper- 
ature increasing  then  decreasing,  the  water  heat  capacity 
increases  so  the  heat  flux  increases  raising  the  wall 
superheat  and   the  heat- transfer   coefficient. 

E.  TURBO -B  TUBE 

Figure  5.  8  shows  a  similar  relationship  as  Figure  5. 5 
but  for  the  Turbo-B  tube.  The  increase  in  wall  superheat, 
for  a  heat  flux  of  40  kW/m2 ,  ranges  from  a  factor  of  1.  04 
for  1  percent  oil  to  a  factor  of  1.  37  for  6  percent  oil.  The 
Turbo-B  tube  remained  active  over  the  full  active  length  of 
the  tube,  as  shown  in  Figure  5.9.  The  apparent  decrease  of 
nucleation  sites,  as  shown  by  the  High  Flux  tube,  did  not 
occur  with  the  Turbo-B  tube.  As  discussed  pertaining  to  the 
High  Flux  tube,  the  increase  in  the  wall  superheat  with 
increasing  oil  concentration  is  probably  due  to  the 
increased  surface  tension  and  viscosity  of  the  mixture 
within  the  cavities.  Also,  Boiling  runs  for  constant  water 
velocity  and  varying  water  temperature  were  performed  on  the 
Turbo-B  tube  and  are  also  shown  in  Figure  5.  10.  The  same 
hysteresis  motions  observed  with  the  High  Flux  tube  are 
observed  with  the  Turbo-B  tube. 

F.  COMPARISON  OF  BOILING  HEAT-TRANSFER  COEFFICIENTS 

The  ratios   shown  in   Figure  5. 11   are  h/h    for  the 
Turbo-B  and  High   Flux  tubes  and  h/hs   for   the  smooth  tube. 
All  of  the  ratios  were  taken  at  a  heat  flux  of  40  kW/m2.  The 
High  Flux  tube  shows   the  most  significant  effect  due  to  oil 
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(a) 


(  b ) 


,5'jre  5.9.  Photographs  of  3 oiling  from  Turbo -3 
Surface:  (a)  0"£,  (bj  Z79,    and  (,  c  j  6'Ji 

Oil  Ma 33  Porcencn. 

Note:  3ar  in  front  of  boiling  tube  is  an 
auxiliary  nearer. 
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concentration  which  is  due  to  the  reasons  previously 
discussed.  The  h/hs  ratio  for  the  High  Flux  tube  decreased 
by  65  percent.  For  the  Turbo-B  tube  the  rate  of  degradation 
of  the  boiling  heat- transfer  coefficient  is  less  pronounced. 
The  Turbo-B  tube  shows  a  decrease  of  24  percent  at  6  mass 
percent  oil.  The  larger  cavity  size  of  the  Turbo-B  tube, 
compared  to  the  High  Flux  tube,  results  in  smaller  oil 
concentrations  in  the  oil  rich  layer.  This  observation  may 
explain  the  smaller  degradation  experienced  by  Turbo-B  tube. 

G.   OVERALL  HEAT-TRANSFER  CHARACTERISTICS 

The  overall  heat-transfer  coefficient  versus  decreasing 
water  velocity  for  oil  mass  percents  of  0,  2,  and  6  for  each 
boiling  tube  is  shown  in  Figure  5.12.  For  pure  R-114  with  a 
water  velocity  of  2  m/s,  the  Turbo-B  tube  outperforms  the 
High  Flux  tube  by  a  factor  of  1. 13  and  the  smooth  tube  by  a 
factor  of  4.  While  the  High  Flux  tube  outperformed  the 
Turbo-B  tube  based  on  the  outside  heat- transfer  coefficient, 
the  reverse  is  true  for  the  overall  heat-transfer  coeffi- 
cient. This  relationship  in  U0  is  due  in  part  to  the 
increased  internal  enhancement  of  the  Turbo-B  tube. 
Comparison  of  the  Sieder-Tate- type  constants  for  the  Turbo-B 
and  High  Flux  tubes  (see  p.  34)  shows  that  the  Turbo-B  tube 
has  a  17  percent  greater  inside  coefficient  than  the  High 
Flux  tube.  Both  of  the  enhanced  tubes  dramatically  outper- 
formed the  smooth  tube.  It  should  be  noted  that  the  increase 
in  pressure  drop  due  to  internal  enhancement  ( compared  to  a 
smooth-interior  case)  was  not  considered  here.  The  pressure 
drop  would  be  an  important  factor  in  finally  selecting  the 
most  economical  tube  type. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

Considering  the  data  gathered  in  this  investigation  for 
boiling  of  R-114  and  oil  mixtures  at  13.8  C  boiling  temper- 
ature, the  following  conclusions  are  reached: 

1.  Based  on  the  overall  heat- transfer  coefficient,  the 
Turbo-B  and  High  Flux  tubes  outperformed  the  smooth 
tube  by  a  factor  greater  than  3,  and  the  Turbo-B  tube 
outperformed  the  High  Flux  tube  by  a  factor  of  1. 13. 

2.  At  a  heat  flux  of  40  kW/m2  with  zero  mass  percent  oil, 
the  outside  coefficients  of  the  High  Flux  and  Turbo-B 
tubes  were  factors  of  14.6  and  6.4  times,  respec- 
tively, compared  to  the  smooth  tube.  These  factors 
decreased  to  7.  0  and  4.  9  for  the  High  Flux  and  Turbo-B 
tubes,  respectively,  with  6  mass  percent  oil. 

3.  The  Turbo-B  and  High  Flux  tubes  showed  Sieder-Tate- 
type  constants  of  0.077  and  0.066,  respectively,  com- 
pared to  the  Sieder-Tate  constant  of  0. 027  for  long 
smooth  tubes.  As  noted  in  the  discussion,  these  values 
may  be  up  to  10  percent  lower  for  internally  enhanced 
long  tubes  than  the  values  for  the  short  tubes. 

4.  While  the  High  Flux  tube  outperformed  the  Turbo-B 
tube,  based  on  the  outside  heat- transfer  coefficient, 
the  High  Flux  tube  is  more  susceptible  to  the  presence 
of  oil. 


B.   RECOMMENDATIONS 

Based  on   the  results  obtained  from   this  investigation, 
the  following  recommendations  are  made: 

1.  The  investigation  of  water-heated  tubes  should  be 
upgraded  to  use  long  tubes.  This  would  decrease  the 
entrance  effects  that  short  tubes  experience  and  it 
would  come  closer  to  duplicating  actual  operating 
conditions.  Also,  the  pressure  drop  effect  previously 
mentioned  in  Chapter  V  (RESULTS  AND  DISCUSSION]  could 
be  investigated  in  the  analysis  of  the  long  tubes. 

2.  Studies  of  the  boiling  heat- transfer  coefficients 
should  be  expanded  to  tube  bundles.  The  interaction  of 
tubes  within  the  tube  bundle  is  a  major  factor  in  the 
analysis  of  heat- transfer  data.  The  boiling  heat- 
transfer  coefficient  may  be  strongly  influenced  by  any 
foaming  action  within  the  bundle. 

3.  Further  investigations  should  be  conducted  on  internal 
tube  enhancements  due  to  the  major  thermal  resistance 
being  the  internal  resistance. 

4.  Modifications  to  the  condensing  side  of  the  present 
apparatus  would  allow  for  studies  of  short  water- 
heated  tubes  at  lower  boiling  temperatures.   The  major 
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modification  would  involve  only  the  increase  of  the 
condenser  drain  piping  to  handle  increases  in  the  flow 
rate  of  the  refrigerant  liquid. 
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APPENDIX  A 
DATA  REDUCTION  PROGRAM 

The  data  reduction  program,  DRP4 ,  used  for  this  investi- 
gation is  listed  below.  A  brief  description  of  the  major 
sections  is  given  in  Table  3,  which  is  followed  by  an  indi- 
vidual line  listing  of  DRP4. 
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TABLE  3 
DRP4  MAJOR  SECTIONS 


Line  Number 

10-80 

80-3690 


3700-4165 


4170' 
4225 

4940 

6180 

6415 

6495- 


4220 
4935 

6175 

6410 

6490 

7280 


7285-7390 
7395-9640 
9645-9725 

9730-9755 
9760-10200 

10205-10360 
10365-10475 


Description 

Selection  options  (i.e. 
data,  plot , etc.  ) 
Sub  Program  Main 
-boiling  tube  dimension 
-select  desired  operati 

conditions  (i.e.  initi 

mode,  flow  rate,  Tsat, 
-convert  emf ' s  to  tempe 
-compute  contact  resist 

(electrically  heated  m 
-compute  various  water 
-various  heat  transfer 

calculations(Q,  LMTD, 
-compute  various  Freon 

properties 
-compute  natural  convec 

transfer  coefficient  f 

unenhanced  ends(iterat 

procedure ) 
-compute  heat  loss  for 

ends 
-compute  actual  heat  fl 

boiling  coefficient 
-record  and  store  data 
Functions  for  curve  fit 
various  Freon  114  prope 
Function  for  polynomial 
Sub  Program  Poly 
-calculate  nth  order  po 
Sub  Program  Plin 
-plot  linear  graphs 
Sub  Program  Stats 
-calculate  standard  dev 
Sub  Program  Coeff 
-plot  cross-plot  files 
Sub  Program  Wilson 
-calculate  modified  Wil 

plot  values  ( see  Appx. 
Functions  for  curve  fit 
various  water  propertie 
Sub  Program  Plot 
-plot  non- linear  graphs 
Sub  Program  Symb 
-print  various  symbols 

plotted  graphs 
Sub  Program  Purge 
-purge  unwanted  files 
Sub  Program  Tdcn 
-calculate  water  temper 

rise  corrections  due  t 

pressure  drop  between 

thermopile  locations 
Sub  Program  U  plot 
-print  O   files 
Sub  Program  Select 
-select  various  sub 

program  options 


process 

s 

ng 

ation 
etc.  ) 
ratures 
ance 

ode  only) 
properties 

U  ,  etc.  ) 
114 

tive  heat- 

or 

ive 

unenhanced 

ux  and 


s  of 
rties 
curve  fit 

lynomial 


iations 


son 

C)f 
s  of 

s 


on 


ature 
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":       '    DP'E:  Qc'toDer    "3.    '33* 

20      !    REVISED:        ftoni    37.    1986 

30      '    COM   /Idp/    Idp 

35  PRINTER    IS    1 

40  CALL  Seiect 

*5  INPUT  "WANT  TO  SELECT  ANOTHER  GPTI0N  ( 1 -Y .0-N>?" . Isei 

50  IF  Isel-1  THEN  GOTO  40. 

55  BEE3 

50  3EE° 

65  PRINTER  IS  1 

70  PRINT  "DATA  COLLECTION/REPROCESSING  COMPLETED" 

75  END 

80  SUB  Main 

35  COM  /lap/  IdD 

90  COM  /Cc/  C<7) .leal 

95  COM  /Wil/  D2.Di .Do .L.Lu.Kcu 

100  DIM  Emf  <!2)  ,T<12>  .Dla<9>  .D2a(9)  .Dial  3)  .Doa<9>  ,La(3>  .Laa<3>  .KcaaO)  ,Et(  19) 

TnSC4)[ 15] 

105  DATA  0.  10036031 .25727.34363.-767345.3295.78025595.31 

110  DATA  -3247486589. 5. 37638E-1 1  ,-2. 66192E-13. 3. 94073E  +  14 

1 15  READ  C(*> 

120    !    DATA    "Smooth"  .  "H  lgn    Flux" , "Thernoexel-E" . "Thermoexel-HE" 

125  DATA    Smooth.  Hign    Fliuc.TurDo-3.Hign    Fiox    Mod  .  TurDO-3    Mod 

130  READ    TnSC) 

135  PRINTER  IS  701 

140  BEE? 

?45  IF  IdD=4  THEM  PRINTER  IS  1 

150  IF  IdP=4  THEM  GOTO  1230 

155  INPUT  "ENTER  MONTH.  DATE  AND  TIME  ( MM : DD :HH : MM : SS) " , DateS 

160  OUTPUT  709:',TD":Date5 

165  OUTPUT  709:"TD" 

170  ENTER  709:DateS 

175  PRINT 

180  PRINT  "  Month,  date  and  time  :":Date5 

135  PRINT 

!90  PRINT  USING  "1  OX  .  ""NOTE :  Drogram  name  :  DRP4 

195  3EED 

20(J  INPUT    "ENTER    DISK    NUMBER".  Dn 

205  PRINT   USING    "16X  .""Disk    niuiiDer      =    "".ZZ":Dn 

210  3EEP 

215  INPUT  "ENTER  INPUT  MODE  < 0  =  3054A . 1 =F ILE) "  .  In 

220  3EEP 

225  INPUT  "SELECT  HEATING  MODE  <0=ELEC:  1  -vJflTER)"  .  Ihm 

230  BEEP 

235  INPUT  "ENTER  THERMOCOUPLE  TYPE  (0=NEW. i -OLD) ".leal 

241)  IF  Im-0  THEN 

2  45  BEEP 

250  INPUT  "GIVE  A  NAME  "OR  THE  RAU  DATA  PILE" . D2_r i ieS 

255  PRINT  USING  "16X.""New  file  name:  "" . 1 4A" ;D2_f l leS 

"?  C  Q  °  :  7  P  i  =  ?  0 

265  CREATE^BDAT    D2_r  i  le'5  .Si  ze  i 

270  ASSIGN   >?Fiie2    TO    D2_;:ie'E 

275    1 

290  !  DUMMY  "ILE  UNTIL  Nrun  KNOWN 

235  D1_f  i  le'5  =  "DUMMY" 

290  CREATE  3DAT  Dl_f i leS.Si zel 

295  ASSIGN  *File1  TO  D1  FileS 
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3  0  0  OUTPUT  if:  i*i  :Date'5 

205  IF  lhm=0  "HEM 

3! 0  3EED 

3,5  INPUT  "ENTER  NUMBER  OF  DEFECTIVE  TCS  <  0=DEFAULT)" .  Idtc 

320  IF  Idtc-0  THEN 

325  Ldtcl-0 

330  Ldtc2=0 

3  35  PRINT  USING  ,,TSX.""Mo  defective  "Cs  exist 

340  END  IF 

345  IF  Idtc  I  THEN 

250  BEEP 

355  INPUT  "ENTER  DEFECTIVE  TC  LOCATION" .Ldtc! 

360  PRINT  USING  "16X.""TC  is  defective  at  iocation  "" .D" :Ldtc1 

365  Ldtc2=0 

370  END  IF 

375  IF  Idtc=2  THEN 

380  BEE3 

335  INPUT  "ENTER  DEFECTIVE  TC  LOCATIONS" .Ldtc ! .Ldtc2 

390  PRINT  USING  "16X.""TC  are  defective  at  locations  ""  ,  D  .  4X  .  D"  ;  Ldtc'  .Ldtc2 

395  END  IF 

400  IF  Idtc>2  THEN 

405  BEE? 

410  PRINTER  IS  1 

4  15  BEEP 

420  PRINT  "INVALID  ENTRY" 

425  PRINTER  IS  701 

430  GOTO  310 

435  END  IF 

440  END  IF 

445  OUTPUT  iFilel :Latc!  .Ldtc2 

450  !  Im=1  option 

455  ELSE 

460  3EE3 

465  INPUT  "GIVE  THE  NAME  OF  THE  EXISTING  DATA  FILE" , D2_f i le5 

470  PRINT  USING  "16X.""0id  file  name:  "" .  1  4A" : D2_f l leS 

4  75  ASSIGN  »Fiie2  TO  D2_i'ileS 
480  ENTER  *File2:Nru.n 

-J85  ENTER  *Fiie2:DoldS 

490  PRINT  USING  "16X.""xhis  data  set  taxen  on  :  "" , ! 4A" : Do idS 

495  ENTER  aFile2:Ldtc; .Ldtc2 

50  0  IF  Ldtc!>0  OR  Ldtc2>0  THEN 

505  PRINT  USING  "  *  6X  .  ""Thernocoup  ies  were  defective  at  locat ions : "" .2< 3D . 4X) 

Ldtc!  ,Ldtc2 

5!0  END  IF 

515  ENTER  aFile2:Itt 

520  END  IF 

525  !  IF  Im-fl  AND  Ihm«!  "HEN  '595 

530  3EE° 

5  35  INPUT  "WANT  TO  CREATE  A  PLOT  FILE?  (0-N. i=Y>" , Iplot 
540  IF  Ipiot=1  THE,' I 

545  3EE3 

550  INPUT  "GIVE  NAME  FOR  PLOT  FlLE".P_f lieS 

555  CREATE  3DAT  P_fiieS.4 

560  ASSIGN  i'PLot  TO  P_fiie5 

565  END  IF 

570  IF  Ihm-1  THEN 

575  BEE3 

580  INPUT  "WANT  TQ  CREATE  iJo  FILE"1  (0-N .  ;  -Y  >".Iuf 

535  IF  Iu.f-1  'r^'i 

5  30  3EE? 
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555  INPUi  "ENTER  jo  rILE  NAME"  ,!Jf  i  ie5 

500  CREArE  3DAT  Jf  :  leS  . - 

b05  ASSIGN  JUfile  TO  UfiieS 

610  END  IF 

6!5  BEE? 

629  INPUT  "WANT  TD  CREATE  Re  FILE''  (  0=N.  !  =Y>"  ,  Ire 

S25  IF  Ire-1  THEN 

S30  BEE? 

555  INPUT  "ENTER  Re  ~ILE  NAME" .Rof ileS 

640  CREATE  BDAT  RefiieS.10 

645  ASSIGN  .-'Refile  TO  RefileS 

S50  END  IF 

555  END  IF 

660  PRINTER  13  1 

665  IF  lm=0  THEN 

670  BEEP 

675  PRINT  USING  "4X. ""Select  tufie  numoer 

680  IF  Ihm-0  THEN 

635  PRINT  USING  "5X.""Q   Smooth  4  men  Rer 

690  PRINT  USING  "SX.""1   Smooth  4  men  Cu  (Press/Slide) 

635  PRINT  USING  "6X.""2   Soft  Solder  4  inch  Cu 

700  PRINT  USING  "6X.""3   Soft  Solder  4  inch  HIGH  FLUX 

705  PRINT  USING  "5X.""4   Uieland  Hard  3  inch 

710  PRINT  USING  "SX.""5   HIGH  FLUX  3  men 

715  PRINT  USING  "6X.""5   GEWA-K  19  Fms/m 

720  PRINT  USING  "5X.""7   GEWA-K  26  Fms/m 

725  PRINT  USING  "SX.""8  GEWfl-T  19  Fms/m 

730  PRINT  USING  "6X.""9   GEWA-T  26  Fins/m 

735  ELSE 

740  PRINT  USING  "SX.""0   Smooth  tube 

745  PRINT  USING  "6X.""1   Hign  Flax 

750  PRINT  USING  "6X.""2   Taroo-3 

755  PRINT  USING  "5X.""3   High  Flax  Mod 

760  PRINT  USING  "6X.""4   Taroo-B  Nod 

765  END  IF 

770  INPUT  I-tt 

775  OUTPUT  AFiiel : It t 

730  END  IF 

735  PRINTER  IS  701 

790  IF  Itt<10  THEN  PRINT  USING  "1  6X  .  ",,TuDe  Number:    "".D":Itt 

795  IF  Itt>9  THEN  PRINT  USING  "16X.""TuOe  NajnDer:    "".DD'Mrt 

300  IF  rhm-1  'HEN  PRINT  USING  "16X .""Tube  Type:      "" . !5A" : TnS< I tt) 

305  BEE? 

310  INPUT  "ENTER  OUTPUT  VERSION  (0-LONG. 1 -SHORT. 2-NONE)" , Iov 

315  BEE? 

320  INPUT  "SELECT  (0=LIQ.  I  -VAP  .2=<LIQ+VAP>/2>". Ilqv 

g  25  ' 

330  !  DIMENSIONS  FQ  TESTED  TUBES 

335  .!  ELECTRIC  HEATED  MODE 

840  !  D1 =D lameter  at  thermocouple  positions 

345  DATA  . 01 1 1125 . .01 1 1125 . . 01 1 1 125 . .0129540 . .012446. .0129540. .0100365 

350  DATA  .0100965.  .01157.  .01 157 

355  READ   D1a<*> 

860  D1=D1a(Itt) 

365    ' 

370  !  D2=Diameter  of  tesr  section  to  the  base  of  fins 

375  DATA  .015375.  . 01 53  75 .. 0 1 53  75 . .015324.  .015875.  .015824.  .01270 

330  DATA  .0127.  .0133. .0128 

335  READ  D2a<*) 

390  ! 

395  !  Di=Inside  diameter  of  anennancea  ends 
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3UU 

905 

9»0  ! 

315  ! 

'320 

.0122 

'325 

920  ! 

335  ! 

340 

945 

'350  ! 

955  ! 

960 

365 

970  ! 

975  ! 

980 

985 

990 

9°5  ' 
000! 
005! 
010 
015 
020 
025 
020 
025 
040 
045 
050 
055 
060 
065 
070 
075 
080 
085 
0  90 
095 
100 
105 
110 
115 
120 
!25 
130 
135 

14Q 

145 
150 
155 

•so 
;tv 

;S5 

170 
'  175 
'130 
1  135 


J  n  :  -I   .  'J  I  ,_  ;  ,  .  u  .  _  ,  ,  .Old  ,  ,JIOl,  .UiL,  ,  .v.  Jt.  i-J 

READ  Dia(*> 

Do=0ut3iae  aiaineTer  of  .inennancea  ends 

DATA  .015375.  .015375.  .015875. .015324. .0153  75 


.0 


01270. .01270. .01331 


0762.  .0762. .076; 


READ  Doa<*> 

L='_engtn  of  enhanced  Burrace 

DATA  .1016. .1016. .1016. .1016. .2032. .2032. .2032 

READ  La^~> 

Lu=Length  oi    unennanced  surface  at  the  ends 
DATA  .0254. .0254. .0254. .0254. .0762. .0762. .0  762 
READ  Lua(*> 

<cu=Thernal  Conductivity  of  tube 

DATA  398  .  344 . 344 . 45 . 344 . 45 ,344 . 844 . 398 . 398 

READ  Kcua<*> 

IF  Ihm=1  THEM 

Data  statements  for  'jater  heating  mode 

DATA  0  .  0 1 5375 . 0  .  0 1 5375 . 0 . 0 1 63 . 0 . 0 1 28 . 0 . 0 1 69 . 0 . 0 . 0 . 0 . 0 

READ  D2a<*> 

DATA  0.0127.0.0127.0.0145.0.0127.0.0145.0.0.0.0.0 

READ  Dia<»> 

DATA  0.015875.0.015875.0.0163.0.015375.0.0169.0.0.0.0.0 

READ  Doa<-> 

DATA  0.3048.0.2043.0.3048.0.3048.0.3048.0.0.0.0.0 

READ  La<*> 

DATA  0  .  0254 . 0 . 0254 . 0 . 0254 . 0 . 0254 . 0 . 0254 .0.0.0.0.0 

READ  Lua(») 

DATA  398.45.398.45.398.0.0.0.0.0 

READ  Kcua<») 

END  IF 

D2=02a(Itt) 

Di=Dia< Itt) 

Do=Doa( Itt) 

L  =  La( Itt) 

Lu=Lua< Itt) 

Kcu=«'.cua(  Itt) 

Xn=.3 

Fr=.3 

IF  Itt-0  THEM  Cf-1 .  70E-9 

IF  1 1 1 > 0  THEN  Cf  =  3  .  7037E  +  1 0  . 

A=PI*<Do  2-Di  '2)/4 

p=pT*nQ 

BEEP 

INPUT  "XUBE  INITIATI0N  10  DE.  <'=H0T  WATER  .2=STEAM.3-CGLD  WATER)"  .  It . 


IF  Itim-1  THEN  PRINT  USING  "16X.""Tube  Initiate 
IF  Itim-2  THEM  PRINT  USING  "1 6X .  ""T.iDe  Initiate 
IF  Itim=:-:  THEN  PRINT  USING  "16X.""Tube  Initiate 


Hot  Water 

Stean ' 

Cold  Water' 


INPUT  "TEMP/vEL  MODE:  < 0  =  7-CQNST . V-QEC ; !  =  "-DEC . 7-CONST  :  2  =  T- INC  .  '/-CONST) " 


IF  Itv-0  THEM  PRINT  USING  " 1 6X  .  ""Temp/Ve i  *oae:  Constant/Decreasing' 

IF  Itv=1  THEN  DRINT  USING  "16X . ""Temp/Vei  looe:  Decreas  i  na/Constan  t ' 

IF  Etv=2  THEM  PRINT  USING  "16X  .""Temp/Vel  floce:  Increasing/Constant' 

INPUT  "WANT  TQ  RUN  WILSON  3U3T?  ( 1 »Y ,0-N>"  .  Iuu  1 
IF  Ihm=!  AND  luui  =  0  THEN 
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•  "30  IF  Itt-0  THEM  1\«.022 

1195  IF  Itt-1  OR  *tr=3  "HEN  :  i  =  .  i'J  5  - 

!200  IF  Itt-2  GR  Ett-a  rHEN  Ci-.062 

1205  BEE? 

1210  INPUT  "ENTER  21  < DEF :  M-.032.HF-.059.TB- .062>".Ci 

1215  PRINT  USING  "1 6X . ""3 ieaer-Tate  

1220  PRINT  USING  "16X.,,M   Constant    =  ,,,,.Z.4D":Ci 

1225  END  IF 

1230  IF  Ihm=!  AND  Em- I  AND  luii-l  THEM 

1225  IF  Itt-0  THEM  Ci-.032 

1240  IF  Itt-1  OR  Itt=3  THEN  Ci=.059 

1245  IF  Itt-2  OR  Itt-4  THEM  Ci=.062 

1250  ASSIGN  WileZ    TQ    * 

1255  CALL  Milson<Cf .Ci) 

1260  ASSIGN  SFiie2  TO  D2_?iie5 

12S5  ENTER  '*Fi le2:Nru.n .DoldS .Ldtd  ,Ldtc2  . 1 1 1 

1270  END  IF 

1275  Nsud=0 

1230  IF  Idp-4  THEN  Ihm=! 

1235  IF  Ihm-1  THEM  Nsub-8 

1290  Ntc=b 

1295  J=i 

1300  Sx-0 

1305  Sy=0 

1310  Sxs=0 

1315  Sxy-0 

1320  Repeat:  ! 

1325  IF  lm=0  THEM 

1330  Dtld=2.22 

1335  Ido  =  2 

1340  ON  <EY  0.15  RECOVER  1320 

1345  PRINTER  IS  1 

1350  PRINT  USING  "4X .""SELECT  OPTION 

1355  PRINT  USING  "SX  .  ""0=TAKE  DATA 

13S0  IF  Ihm-Q  THEM  PRINT  USING  "SX.,#,,1=SET  HEAT  ^LUX 

13G5  IF  Ihm-1  THEN  PRINT  USING  "SX.""1=SET  WATER  PLOW  RATE' 

1370  PRINT  USING  "GX.""2=SET  Tsat 

1375  PRINT  USING  "4X.""N0TE:  KEY  0  -  ESCAPE 

1330  BEE? 

1385  INPUT  Ido 

1390  IF  Ido>2  THEM  Ido  =  2 

1395  IF  ldo=0  THEN  2145 

1400! 

1^05!  LOOP  TO  SET  HEAT  FLUX  OR  FLOWMETER  SETTING 

1410  IF  Ido-1  THE.M 

Ml  5  IF  Ihm-0  "HEM 

1420  OUTPUT  709:"AR  AF12  AL13  VR5" 

M25  BEE3 

1430  INPUT  "ENTER  DESIRED  0<jd"  .Docd 

1435  PRINT  USING  "4X . ""DESIRED  Qdo   ACTUAL  QdD 

1440  Err-1000 

M^5  FOR  1  =  1  TO  2 

1*50  OUTPUT  709: "AS  SA" 

1455  Sum=0 

1460  FOR  J 1  =  '  TO  5 

1465  ENTER  709:E 

M70  Sum-Sum+t 

14  75  NEXT  Ji 

1430  IF    r-!    THEN   Voit-Sum*5 

1485  IF    1=2    THEN   AmD=E 
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■  .3  an 

next  : 

!495 

Aaap=*^oi  t*rtriD/  (PI*D2*U' 

1500 

IF  ABS<Aqap-Daap>>Err  THE! 

1 

1505 

IF  Aaap>QaaD  "HEM 

1510 

BEED  4000  .  .2 

1515 

3EE?  4000.  .2 

1520 

BEE°  4000.  .2 

1525 

ELSE 

1530 

BEE3  250. .2 

!535 

BEE?  250.  .2 

1540 

BEEP  250, .2 

15^5 

END  IF 

1550 

PRINT  USING  "4X.MZ.3DE.2X, 

NZ.2DE" 

:  Dadp .  Aqd( 

1555 

WAIT  2 

1560 

GOTO  1445 

1565 

ELSE 

1570 

BEE? 

1575 

PRINT  USING  "4X.MZ.3DE.2X, 

.MZ.3DE" 

1  :  DadD  .  Aadi 

1530 

Err=500 

1535 

WAIT  2 

1530 

GOTO  1445 

1595 

END  IF 

1600 

ELSE 

1605 

BEE? 

1610 

INPUT  "ENTER  FL0UMETER  SEl 

"TING".Fm< 

16*5 

GOTO  1350 

1620 

END  IF 

1625 

END  IF 

1630' 

1535! 

LOOP  TO  SET  Tsat 

164Q 

IF  Ido  =  2  THEN 

'645 

IF  Ikdt-1  THEN  1670 

1650 

BEE3 

1655 

INPUT  '"ENTER  DESIRED  TSat' 

'.Dtid 

1660' 

PRINT  USING  "4X.""  DTsat 

ATsat 

Rare 

1655 

Ikdt-1 

1670 

a  1  d 1  -  0 

1675 

Qld2*0 

1630 

Mn-1 

'635 

Nrs=Nn  *OD  15 

1630 

Nn=Nn+i 

'695 

IF  Nrs-i  THEN 

:  700  _ 

PRINT  USING  ••4X.,•••   Tsat 

Tldl 

Tld2 

T 

out 

1705 

END  IF 

1710 

IF  Ihm-0  THEN  OUTPUT  703:' 

•AR  AF3 

AL'I  VR5" 

1715 

IF  Ihm-i  THEN  OUTPUT  709:' 

"AR  AFO 

AL5  VR5" 

1720 

FOR  I-l  T0  6 

1725 

IF  lhm=0  AND  I>4  THEN  1860 

1730 

Subi*0 

1735 

OUTPUT  709: "AS  SA" 

1740 

FOR  Ji-1  TO  20 

'  745 

ENTER  709:Eiig 

1750 

Sum=Siun+E  i  iq 

!  755 

NEXT  Ji 

1760 

Eliq=Sum/20 

1765 

rid-FNTvsv(El La) 

'770 

IF  1  =  1  THEN  TldWld 

•7^5 

IF  1=2  THEN  TId2=T;d 

■  730 

IF  1=3  THEN  Tv=~ld 

Rat. 


Tslutid         Tinlet  Tpile 
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.  -  ac 

IF  I  =  J  "HE: J  TsumD  =  7id 

■  790 

IF  1=5  THEN  Tinier-Tic 

:  735 

IF  I"-.'  i  HEM  "roixt  =  TId 

1800 

NEXT  I 

1805 

IF  Ihm-I  THEM 

1310 

OUTPUT  709:"AR  AF20  AL20  VR5" 

1315 

OUTPUT  709: "AS  SA" 

1320 

Sum«0 

1325 

FOR  Klc-I  TO  20 

1330 

ENTER  709 :E 

1335 

Sum-Sum+E 

1840 

NEXT  Kk 

1345 

Emf  <  7>=4BS<Sum/20> 

1350 

TDile=Emt (7)/3.96E-4 

1855 

EMD  IF 

1360 

Atld=< Tidl+Tld2)*.5 

1865 

IF  ABS(Atld-Dtid>>.2  THEM 

1370 

IF  Atid>Dtld  THEM 

1375 

BEEP  4000.  .2 

1330 

BEEP  4000.  .2 

1835 

BEE3  4000. .2 

1890 

ELSE 

1895 

BEEP  250.  .2 

1900 

BEE?  250.  .2 

1905 

BEE0  250.  .2 

1910 

EMD  IF 

1915 

Err1-Atld-Old1 

1920 

Oldl-rttld 

1925 

Err2»Tv-Qld2 

1930 

Qld2»Tv 

1935 

IF  Ihm-0  THEM  PRINT  USIMG  "AX. 

,5<MDD. 

DD.2X)":Dtid. 

Tldl  .  Tld2  .  Tv  ,  Tsixrnp 

1940 

IF  Ihm-1  AND  ldp=0  THEN  PRINT 

USING 

"4X.7(MDD.DD, 

2X)";Qtid.Tld1 .Tid2 

ujnp  .  T 

iniet . Tpi le 

1945* 

IF  Ihm=!  AND  Idp=^  THEM  PRINT 

USING 

"4X.5(MDD.DD. 

2X>,3<M3D.DD.2X)":D 

dl  .Tld2.7v.Tsixmp.7in 

1950 

^!AIT  2 

1955 

GOTO  1635 

1360 

ELSE 

1965 

IF  ABS(Atld-Dtid)>. !  THEN 

1970 

IF  Atld>Dtid  THEN 

1975 

BEE?  2000.  .2 

1930 

BEE?  3000. .2 

1985 

ELSE 

1990 

BEE?  800.  .2 

1995 

BEEP J300. .2 

2000 

END  IF 

2005 

Err ' =At id-Q Id ! 

2010 

0  Id  1 »At id 

2015 

Err2- i v-Qld2 

2020 

01d2="v 

2025 

IF  Ihm-0  "HEM  PRINT  USING  "4X 

,5<MDD 

,DD.2X>":Dtid, 

,  T  id  1  ,  T lo2 . Tv . "sumo 

2030- 

■  IF  Ihm=1  THEN  PRINT  USING  "4X 

.5<«DD 

.DD.2X) ,3<M3D 

,  3D. !X)":Dtid. Tld1  . f 

.  T  sixmo  .Tiniet.Tpile. 

2035 

WAIT  2 

2040 

GOTO  1635 

2045 

ELSE 

2050 

BEEP 

2055 

Err1 =Atid-Old1 

2060 

01d1=Atld 

2065 

Err2-Tv-01d2 

Ld.Tl 


[dZ 
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20  70  Qid2*Tv 

2075  IF  Ihm-0  'MEM  PRINT  USING  "4X .5<MDD .DO .2X>" :Dtld . Tid! . Tld2 . Tv. Tsumo 

2080  IF  Ihm-1  THEN  PRINT  USING  "4X  .  3  <  MDD  .  DD  .  2X  )  "  :  D  t  :o  .  T  id  1  .  Tld2  .  Tv/ .  Tsmnc  .  7  in  le  t 

, Td i le . Tout 

2035  WAIT  2 

2090  GOTO  1685 

2095  END  IF 

2 100  END  IF 

2105  END  IF 

2110!  ERROR  TRAP  FOR  Ido  GUT  OF  3QUMDS 

2115  IF  Ido>2  THEN 

2120  3EEP 

2125  GOTO  1350 

2130  END  IF 

2135  * 

2140!  TAKE  DATA  IF  Im-0  LOOP 

2145  IF  Ikol-1  THEN  2165 

2 '50  3EEP 

2155  INPUT  "ENTER  SULK  OIL  '/.".Boo 

2160  Ikoi-1 

2165  IF    Ihm-0    THEN   OUTPUT    709:"AR   AF0   ALU    VR5" 

2170  IF    Ihm-1    THEN   OUTPUT    709:,,AR   AF0    AL5    VR5" 

2175  IF    Ihm-0    THEN   Ntc-12 

2130  FOR    1=1    TO   Ntc 

2185  OUTPUT    709: "AS   SA" 

2190  Sum-0 

2195  FOR    Ji-1    TO   20 

2200  ENTER    709:E 

2205  Smn-Sum+E 

2210  IF    I-<9-Nsu£>)    OR    I  =  (!0-Nsud)    THEN   Et<Ji-l)-E 

2215  NEXT   Ji 

2220  Kdi=0 

2225  IF    I-<9-Nsub>    OR    I=(10-NsubJ    THEN 

2230  Eave=Sum/20 

2235  Sum-0. 

2240  FOR  Jlr  =0  TO  19 

2245  IF  ABS(Et<Jk)-Eave><5.0E-6  THEN 

2250  SLun  =  SLUn--Et(Jlc  ) 

2255  ELSE 

2250  Kdl-Kdl+1 

22S5  END  IF 

2270  NEXT  Jk 

2275  IF  I  =  (9-Nsud)  OR  I-<10-Nsub>  THEN  PRINT  USING  "4X."",<dl  =  "",DD" :Kdi 

2230  IF  Kdl>i 0  THEN 

2235  3EE3 

2290  BEEP 

2295  PRINT  USING  "AX. ""Too  Tuch  scattering  in  data  -  reoeat  data  set 

2300  GOTO  1345 

2305  END  IF 

23i Q  END  IF 

2315  Emt(I)=Sum/<20-Kdl> 

2320  NEXT  I 

2225  IF  Ihm-1  THEM 

2330  OUTPUT  709:"AR  AF20  AL20  VRS" 

2335  OUTPUT  709: "AS  SA" 

2340  Swn=0 

2  3^5  FOR  Kle-|  TO  20 

2250  ENTER  709:2 

2255  S  i.im  =  S  am  +  E 

2360  NEXT  Kk 

2265  Emf<7)=ABS<Sum)/20 
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1:270  END  IF 

2375  IF  Ihm-0  "HEN 

2380  DU"?!Jr  709:,,AR  ^F:2  AL:3  VR5" 

2335  "OR  1-1  TO  2 

2330  OUTPUT  709 : "q^  SA" 

2395  Sun-0 

2400  FOR  Ji-I  rQ  2 

2405  ENTER  709  :E 

2410  Sum-Sum+t 

2415  NEXT  Ji 

2420  IF  1='  THEN  Vr-Sum/2 

2425  IF  1  =  2  THEM  Ir-Sujn/2 

2*30  NEXT  I 

2435  END  IF 

2440  ELSE 

24/i5  IF  Ihm=0-THEM  ENTER  W:  ie2  :  Boo  .  To  IdS  .Emt<  *>  .Vr  .  Ir 

2450  IF  Ihm-1  THEN  ENTER  i>F i ie2 : Bop . To idS . Zmi (*). Fms 

2455  END  IF 

2460! 

2465!  CONVERT  emf*  S  TO  TEMP .VOLT .CURRENT 

2470  Tua=0 

2475  FOR  I=i  TO  Ntc 

2480  IF  ldtc>0  THEM 

2485  IF  I=Ldtc!  OR  I=Ldtc2  THEN 

2490  T<I>=-99.99 

2495  GOTO  2545 

2500  END  IF 

2505  END  IF 

2510  IF  Itt<4  AND  Ehn-0  THEN 

2515  IF  I>4  AND  I<9  THEN 

2520  T<I)=-99.99 

2525  GOTO  2545 

2530  END  IF 

2535  END  IF 

2540  T(I)=FMTvsv(Emt (I)) 

2545  NEXT  I 

2550  IF  Itt<4  THEN 

2555  FOR  1=1  TO  4 

2560  IF  I=Ldtc1  OR  I=Ldtc2  rHEN 

2565  Twa=  i ua 

2570  ELSE 

2575  Tua=Tua*-T(  I) 

2580  END  IF 

2535  NEXT  I 

2590  Tua=Tua/<4-Idtc> 

2595  ELSE 

2600  IF  Ihm=!  THEN  2650 

2605  FOR  1=1  tq  g 

2610  IF    I-Ldtcl    OR    I=Ldtc2    THEN 

2515  Tua=Tya 

2620  ELSE 

2625  Tua=Tua+T< I > 

2530  END    IF 

2635  NEXT    I 

2640  Tdja=Tijia/<8-Idtc) 

2645  END  IF 

2650  Tld=T(9-Nsub) 

2555  Tld2=T<10-Nsub) 

2660  Tlda=< TLd+Tid2; * . 5 

2685  Tv=T( 11 -Nsub' 
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lor "j  :-   :tt<;  -no  ihm-o  :hem 

2675  Tid2--9S.99 

2680  ">/  =  '■('  n  )  +  "<;:))  ,'2 

2635  END  IF 

2630  IF  Ihm-0  THEN  2710 

2695  T3uLfnD  =  7  <  12-Nslld) 

2700  Tmiet  =  T(  !3-Ms>xd> 

2705  Tout  =  r<  1  4-Nsub  > 

2710  IF    Ihm-0    THEN 

2715  flmo-flBS(Ir) 

2720  Voit-ABS<Vr)~25 

2725  0=Voit*Amp 

2730  END  IF 

2735  IF  Itt-Q  AND  Ihm-0  THEM 

2740  Kcu-FMKca(Tua) 

2745  ELSE 

2750  Kcu-<ciia(Itt> 

2755  END    IF 

2760' 

27G5!  -OURIER   CONDUCTION   EQUATION   WITH   CONTACT   RESISTANCE    NEGLECTED 

2770  IF    lhm=0    THEN    Tu=Twa-Q*L0G<D2/D1 > /<2*PI»Kcu*L) 

2775  IF    Ilqv-0    THEN   Tsat=Tlda 

2780  IF    Ilqv=1    THEN    Tsat=< Tida+Tv)» .5 

2735  IF    Hqv-2   THEN    Tsat-Tv 

2790  IF    Ihm-1    THEN 

2795  Tavg  ■■.  7  i.  n  1  e  t 

2300  Grad-37.9853+.104388*Tavg 

2305  Tdroo=ABS(Emt (7) )*! .E+6/< 10»Grad) 

2810  Tavgc-Tinlet-Tdrop* .5 

2915  IF  ABSCTavg-Tavgc>>.01  THEN 

2320  Tavo= ( Tava+Tavqc >* . 5 

2325  GOTO  2900 ' 

2330  END  IF 

2335 ! 

2340!  COMPUTE  WATER  PROPERTIES 

2845  Kui«FNKw<Tavg) 

2350  Maua=»-''IMij.ut  Tavg) 

2955  Cpu-FNCpu<Tava) 

2960  Pru  =  -MP-'j(  Tavq) 

2865  Rhow  =  «rNRhou(Tav-g) 

2970  TUi=Tavg 

2375! 

2380!  Compute  MOOT 

2385  ridot=3.9657E-3+Fms*<3.61955E-3-Fms*<3.82006c-6-Fms*< 1 .23638E-7-Fms*4 . 31397 

E-i0)  )  ) 

2990!  *!dot-Mdot*<  I  .0365-Tinlet*<  1  .36644E-3-Tinlet*5  .252E-6)  )/1  .0037 

2995  K  d  t - 0 

2900  fl=i'1dot»Cpiu*Tdrop 

2905  Lmtd=Tdrop/LQG( <  TiRlet-Tsat)/< Tiniet-Tdrop-rsat ) > 

2910  Uo-Q/<PI*Oo*L*Lmtd) 

29  1 5  Ru  =  Do*L0G <  Do/ D i ) / ( 2 . *<cu) 

2920  Tu-Tsat+Fr*Lmtd 

2925  Vw-!1dot/<Rhou*PI*Oi  '2/4) 

2930!  IF  Kdt=0  THEN 

2935!  Kdt-1 

2940!  rdrop"Tdrop-.004*Vu"2 

2945!  GOTO  3710 

2950!  END  IF 

2955  Rew=Rhow*<Vu*Qi/MuAda 

2960  Hi=Ci*'<'iij/Di*Reuj    .8*pru    (t/3.  )«»<Muuia/FNMimi<  Tun)  >"  .  14 

2965  Tuic=Tavg-0/(PI»Qo,*L*Hi  ' 
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2970  IF  «BS<7ui-tuiC)>. j|  THEN 

.-3  7  5  i  u  i =  <  :  'ij  i  *  •  u  i  c  >  ■*  .  ~j 

2980  GOTO  2960 

2935  END  IF 

2390  Twi  =  < Twi  +  Tuic)*  .5 

2995  Ho=1/<1/Uo-Do/<Di-Hi>-Rw> 

3000  ThetaD=Q/<Ho*?I*Qo*L) 

3005  Tw= Tsar +■ The taD 

3010  ThetaD'Tu-Tsat 

3015  IF  Thetab<0  THEM 

3020  BEE? 

3025  INPUT  "TUALl^TSAT  (0-CONTINUE.  f"END>",Iev 

3030  IF  Iev-0  THEN  GOTO  1325 

3035  IF  lev- I  THEN  3535 

30*0  END  IF 

3045  END  IF 

3050! 

3055!  COMPUTE  VARIOUS  PROPERTIES 

3060  Tfilm=-(Tu)+Tsat>*.5 

3065  Rho  =  »rNRho<Tf  ilm) 

3070  Mu=FNMu(Tf llm) 

3075  K=FNK(Tfilm) 

3030  Cp=FNCp(Tf ilm) 

3085  3eta=>-NBeta<Tf  ilm) 

3090  Hfg=FNHTa<Tsat) 

3095  Ni=Mu/Rho 

3100  AlDha=K/(Rho»Cp) 

3105  Pr-Ni/ftlpha 

3110  Psat'FNPsat(Tsat) 

3115! 

3120!  COMPUTE  NATURflL-CGNVECTIVE  HEAT-TRANSFER  COEFFICIENT 

3125!  FOR  UNENHANCED  END<S) 

3130  Hbar=!90 

3135  Fe-<Hbar»P/(Kcu*A> )  "  ,5*Lu 

31 40  Tanh*FNTanh<Fe> 

31^5  Theta=ThetaD*Tann/ Fe 

3150  Xx  =  <  9  .3  1  *6eta*ThetaD*0o    3**Tann/  (  Fe*N  i  *Alpha>  >'  .'65667 

3155  Yy  =  (U<  ,559/Pr)    <  9/ 1  6 )  > "  <  8/27  ) 

3160  Hbarc=</Do~< .G+.387*Xx/Yy) '2 

3165  IF  ABS<(Hbar-Hbarc)/HbarO>  .001  THEN 

3170  Hbar=(Hbar+Hbarc>*.5 

3175  GOTO  3135 

3130  END  IF 

3135! 

3190'  COMPUTE  HEAT  LOSS  RATE  THR0UG_H  UNENHANCED  END(S) 

3195  Ql=<Hbar-P*Kcu*tt)  .5*Thetao- i anh 

3200  Qc=Q-2*Oi 

3205  As-PI*D2*L 

2210! 

32*5!  COMPUTE  ACTUAL  HEAT  FLUX  AND  BOILING  COEFFICIENT 

3220  Qdp=Qc/As 

3225  HtuDe=OdD/ The tao 

3220  Csf  =  <Cp*ThetaD/'Hfg)/(Odp/<Ma*Hfg)*<  .014/(9. 31»Rho)  ,5>"n/3.)*Pr  1.7) 

?  -->  O  C  I 

3240!  RECORD  TIME  OF  DATA  TAKING 

32^5  IF  Im-0  THEN 

3250  OUTPUT  709:,,TD" 

3255  ENTER  709:ToldS 

32S0  END  IF 

3265  ! 
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run  No 

0  i  i  7. 

Tsat 

Htube 

Qdo 

Thetao 

FMS 

n  r  i  •' 

!J  1 L_  A 

rSAT 

HTUBE 

QDP 

THETAB 

t-270!  JUTPUT    DATA    7Q    ^RIM""5 

327C,  PRINTER    13    70"! 

3230  IF    Iov-0    THEN 

3285  DRIN7 

3290  dRINT  USING  '"OX. ""Data  Set  Number  -  ""  .DDD  .2X  .""Bulk  Oil  ■'  -  ""  .DD  .  D  .5X .  \ 

4A"  ;J.Bop  .  i  oi.JS 

3295  IF  Ihm-0  'HEM 

3300  PRINT  USING  "10X.""TC  No:     !      2      3      4      5      S      7 

1.305"  DRINT  USING  "10X.""7emp  :  ""  ,3<  1 X  .MDD  .  DD>"  :  T(  1 )  ,  T<2>  .  T<3>  ,  T<4>  .  T(5)  .  T(S)  .  T< 
7  )  ,  T  (  3  > 

3310  PRINT  USING  "1QX.""  Tuia.  Tli-qd    TIiaa2   Tvapr   Psat    Tsumo 

3315  PRINT  USING  "1  OX  .2(MDD  .  3D  .  1  X)  .  1  X  .MDD  .  DD  .  1  X  .2(  1  X  .  MDD  .  DD)  .2X  .  MOD  .  D"  :  T'ua.  Tid  . 
Tld2 . Tv . Psat . "sump 

3320  PRINT  USING  "10X.""  Thetab   Htube      Qdo 

3325  PRINT  USING  '"!  OX  ,MDD  .  3D  .  1  X  .MZ  .3DE  .  1  X  .MZ  .3DE"  :  Thetao  .Htuoe  .  Qap 

3330  ELSE 

3335  PRINT  USING  "10X.""  Fms    ^u    Tsat  T  ml    Tdrop  Thetab     q         Uo 

Ho 

3340*  PRINT  USING  "1 OX . 4(2D . DD . ! X > . Z . 3D . : X . DD . DD . ! X . 3<MZ . 3DE . ! X)" :Fms . Vu . "sat . Ti 

nlet. Tdrop. Thetao. 3d 

3345  END  IF 

3350  END  IF 

3355  IF  Iov-t  THEN 

3350  IF  J-l  THEM 

3365  PRINT 

3370  IF  Ihm-0  THEM 

3375  PRINT  USING  "10X. 

3330  ELSE 

3385  PRINT  USING  "10X. 

3330  END  IF 

3395  END  IF 

3*00  IF  Ihm-0  THEM 

3405  PRINT  USING  'M2X  .3D  .4X.DD  .2X  .MDD.DD  .3nX.MZ.3DE>" ;  J  .Boo.  Tsat  .Htube.Qdp  .  The 

tab 

3410  ELSE 

3415  PRINT  USING  "12X.3D.4X  .DD.2X  .MDD.DD.3C  IX  .MZ.SDE)"  :Fms  .3op  . Tsat. H tune. Qdp.T 

hetao 

3420  END  IF 

3425  END  ~F 

3430  IF  Im-0  THEN 

3^35  BEED 

3440  INPUT  "QK  TO  STORE  THIS  DATA  SET  C T-Y ,0-N>?".Qk 

3445  END  IF 

3450  IF  Qk-1  3R  Im- I  THEN  J-J-H 

3455  IF  Qk-1  AND  Im-0  "HEN 

3-60  IF  Ihm-0  THEN  OUTPUT  ?FileI :3op.ToidS.Emf (*) ,Vr .Ir 

3465  IF  Ihm=1  THEM  OUTPUT  <W i ie1 ;Bop . ToidS ,Emf < *) .Fms 

3470  END  IF 

3475  IF  Iu.f-1  THEM  OUTPUT  Sljf  l  le  ;Vu».Uo 

3430  IF  Ire-1  THEM  OUTPUT  4Ref i le :Fns .Rew 

3-85  IF  (Im- I  OR  Qk-1)  AND  loiot-i  "HEN  OUTPUT  W lot :Qdp. Thetab 

3430  IF  Im-0  THEM 

3^95  SEE3 

3500  INPUT  "VIILL  THERE  3E  ANOTHER  RUN  (  !  -Y.  0-N>?"  .Go_on 

3505  Nrun=J 

35 !0  IF  Go_on=0  THEN  3535 

35' 5  IF  Go_on<>0  THEN  Repeat 

3520  ELSE 

3525  IF    J<Mrun+l    "HEN    R<-DPat 

75  30  END    :r 
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?.^35  It    i<n-.j  THEN 

3540  3EE- 

?.5a5  PRINT  USING  "I  OX. ""'MOTE:  '"".ZZ.""  data  runs  were  scored  m  file  ""  .  :  0A"  :  J- 

1 .D2_f ileS 

3550  ASSIGN  sFilel  TO  * 

3555  OUTPUT  ■SFile2:Nran-i 

3560  ASSIGN  I'Fi  lei  TO  Dl_fiieS 

3565  ENTER  'IFi  lei  :DateS.Ldtc1  .Ldtc2.  Itt 

3570  OUTPUT  *Fi  ie2  :  QateS  .Ldtc  i  .Ldtc2  .  Itt 

3575  FOR  1-1  TO  Nrun-1 

3530  IF  Ihm-0  THEN 

3585  ENTER  SFilel ;3op.ToldS.Emf (») .Vr . Ir 

3590  OUTPUT  $FUe2:Bop  .ToldS.Emf  <*>  .Vr  ,  Ir 

35S5  ELSE 

3G00  ENTER  »Fi lei :Bop  .  ToldS .Emf <*> ,Fms 

3605  OUTPUT  $Fi  le2 :Bop . ToidS .Emf < *> .Fms 

3610  END  IF 

3615  NEXT  I 

3620  ASSIGN  »Filei  TO  * 

3625  PURGE  "DUMMY" 

3630  END  IF 

3635  BEEP 

3640  PRINT 

3645*  IF  Iplot-1  THEM  PRINT  USING  "TOX .""NOTE:  "".ZZ.""  X-Y  pairs  were  stored  m 

p  iot  data  file  "" , 1 

3650  ASSIGN  :«F;le2  TO  * 

3655  ASSIGN  *Plot  TO  * 

3660  IF  Iuf-1  THEN  ASSIGN  ?Ufile  TO  * 

3665  IF  Ire=1  THEN  ASSIGN  SRefile  TO  * 

3670  CALL  Stars 

3675  BEE3 

2630  INPUT  "LIKE  TO  PLOT  DATA  ( 1 = Y , a-N) 7" , Gk 

3635  IF  Gk  =  1  THEN  CALL  Plot 

3630  SUBEND 
3695' 

3700!  CURVE  FITS  OF  PROPERTY  FUNCTIONS 

3705  DEF  FNKca(T) 

3710!  0FHC  COPPER  25  0  T0  300  < 

3715  Tk  =  ^273.  15     !C  TO  K 

3720  <  =  <j34-.  M2*fk 

3725  RETURN  K 

3720  FMEMD 

3735  DEF  FNMu(T) 

37^0!  170  TO  360  K  CURVE  FIT  OF  VISCGUSITY 

3745  Tk=T+273.15      !C  TO  K 

3750  Mu=EXP< -4. 4636+ ( 1 011 ,47/Tk) )*! .OE-3 

3755  RETURN  Mu 

3760  FMEMD 

3765  DEF  FNCp<T) 

3770!  130  TO  400  '<  CURVE  FIT  OF  Cp 

3775  Tk=T+27o. 15     \r    -p  < 

3730  Cp-.  40 188  +  1  .5500  72-  3*~k  +  i  .5  1  494E-5-T><  2-6 .S7353E-19*Tk  '3 

3795  Cp-Cd*1Q0Q 

3790  RE 'URN  Cp 

3795  FNEMD 

330  0  DEF  -NRho<T) 

3805  Tk=T+273. 15     !C  to  < 

3310  X=!-H  .3*Tk/753.'35)   !K  "0  R 

3315  Ro  =  36. 22+6!  .14641 4*;<"<  T  /3>  +  16  . 4 T 30 1 5*X-*- 1  7  ,476338»X  \5+t  .1  !9828*X  "2 

3320  Ro=Ro/ . 062423 

. 5  RETURN  Ro 


Ik 


2320 
3005 

3940 

3345 

3350 

3955 

3'3S0! 

3365 

3370 

3875 

3380 

3335 

3890 

3395 

3900 

3905 

3910 

3915 

3920 

3925 

3930 

3935 

3940 

3945 

2950 

3955 

3960 

2965 

3970 

3975 

3980 

3985 

3990 

3295 

4000 

4005 

4010 

4015 

4020 

4025 

4030! 

4025 

4040 

4045 

4050 

4055 

4060 

4065 

4070 

4075 

4030 

4035 

4090 

4095 

4100 

4  105 

41  10 

4  !  15 

■4  120 

4125 


) 


leal 


fnend 

DEr  FNPr<T) 

Pr-rMCo<T>*FMMa<  T>/FNK< 

RETURN  °r 

"MEND 

DEF  FNK(T) 

T<3S0  K  WITH  T  IN  C 

K-  .071 -.000261 *T 

RETURN  K 

FNEMD 

DEF  FNTanMX) 

P=EXP<X) 

Q-t/P  ' 

Tanh-<P-Q>/<P+Q> 

RETURN  Tann 

FNEND 

DEF  FNTvsv(V) 

COM  /Cc/  C(7) 

T-C(O) 

FOR  I=i  TO  7 

T»T+Cd)»<r  I 

NEXT  I 

IF  Icai=1  THEN 

T=T-6.7422934E-2+T»<9.0277043E-3-T-<-9.32599l7E-5)) 

ELSE 

T-T+8.62G897E-2+T»<3.76I99E-3-T*5.0639259E-5) 

END  IF 

RETURN  T 

FNEND 

DEF  FNBeta(T) 

RoD=frNRho(  r  +  .1> 

Rom=FNRho(T-. 1 ) 

Beta=-2/(Rop-,-Rom>*<RoD-Rom)/.2 

RETURN  Seta 

FNEND 

DEF  FMHfg(T) 

Hfg-1 .3741344E+5-T*<3.309436iE-<-2+T»l ,2!65I43> 

RETURN  Hfg 

FNEND 

DEF  FMPsat(Tc) 

0  TO  30  aea  F 

Tf-1  .3-Tc+32 

Pa=5.345525^T 

Pq=Pa-;4.7 

IF   Pg>0    THEN 
Psat=p3 
ELSE 

P-3at=°'3-29.'32/i4 
END  IF 
^E'URM  3sat 
FNEND 

DE~  FNHsmooth<X.3op 
DIM  A<5) .3(5) .C(5) 
DATA  .20526. .25322 
DATA  .74515. .72292 
DATA  .41032.  .  1  7725 
DATA  .71403. .72913 
READ  A(*).8<*>.C(* 
IF  Boo<S  THEN  I=3op 
IF  3od=6  THEN  1=4 


CURVE 

FIT  OF 

Psat 

*<  . '5352082+' 

"f*< ; 

.43<J 

0963E- 

■3+ 

rf 

*9 

5  1 

506 

1 

*=^SIG 

-=m 

Hg 

1E-6)) 


3 . Isa" ) 

.  D  <  5 ) 

.  .319048 

. .55322 

.79909 

' .00253 

.  .73189. 

.71225. 

.634  72. 

^-,?7 

. .25142. 

.54306. 

31916. 

.0845 

.  .72565. 

.696691 

. .66536; 

".  .51339 

> .Q<*) 
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-:'T'J  '.-    3op=:0  THEN  1=5 

41  35  IF  :;ar  =  '  "HEN 

4!  4ij  Hs=EXP<A<  I>+3<  I>*'_OG(X)  ) 

4145  ELSE 

4  150  Hs=£XP(C<I)+D< I)*LQG(X) ) 

4155  END  IF 

4  160  RETURN  Hs 

4165  FNEND 


5495 

SUB  Uilson(Cf .Ci > 

6500 

COM  /Mil/  D2.Di.Do.L.Lu.Kcu 

6505 

DIM  Emt< 12) 

5510! 

WLISON  PLOT  SUBROUTINE  DETERMINE  CF  AND  CI 

6515 

BEE? 

6520 

INPUT  "ENTER  DATA  FILE  NAME".FileS 

6525 

BEED 

6530 

PRINTER  IS  1 

6535 

PRINT  USING  "4X. ""Select  ootion: 

6540 

PRINT  USING  "4X.""   0  Vary  Cf  ana   Ci 

55^5 

PRINT  USING  "4X.""   1  Fix   Cf  Vary  Cl 

6550 

PRINT  USING  "4X.""   2  Vary  Cf  Fix   Ci 

655! 

PRINT  USING  "4X.""   3  Fix   Cf  Fix   Ci 

6555 

INPUT  "ENTER  OPTION" . Icf ix 

5560 

PRINTER  IS  701 

6565 

IF  Icfix-0  THEN  6535 

6570 

IF  lcfix>0  THEN  3EEP 

6575 

IF  Icfix-1  THEM  INPUT  "ENTER  Cf'.Csf 

6530 

IF  Icfix=2  THEM  INPUT  "ENTER  CF'.Ci 

5531 

IF  Icfix-3  THEN  INPUT  "ENTER  Cf.  Ci".Csf.Ci 

5535 

PRINTER  IS  1 

5590 

INPUT  "Want  To  Varv  Coef  f  ?<  1  =Y  .  0  =  N  )  "  .  Iccoet 

6595 

IF  Iccoet -1  THEN  INPUT  "ENTER  C0EFF".R 

6600 

PRINTER  IS  701 

5605 

IF  Icfix-0  OR  Icfix-2  "HEN  Cfa=.004 

6610 

IF  Icfix-1  THEN  CFa-Csf 

5515 

Cia=Ci 

5620 

Xn=  .3 

5625 

Fr-.3 

6630 

Jj  =  0 

5535 

Rr=3. 

5640 

IF  Iccoet ='  THEN  Rr  =  R 

6545 

PRINTER  IS  1 

6650 

PRINT  Do.Di  .;<cu 

6655 

ASSIGN  *File  TO  FileS 

6660 

ENTER  ?Fi le:Nran. DateS. Ldtc! .Ldtc2.Itt 

66o5 

Ru=Do-lQG<Do/Di )/(2*Kca> 

5570 

Sx-0 

6675 

Sy  =  0 

5630 

Sxy  =  0 

5635 

Sx2-0 

5590 

Sy2  =  0 

6695 

FOR  1=1  TO  Nrun 

5700 

ENTER  "File: 3od . To IaS . Emf (»> . Fms 

6705! 

CONVERT  EMF  Z    70  "EMPERATURE 

5710 

FOR  J-1  TO  5 

6715 

T<J)=FNTvsv(Emf (J)  ) 

6720 

NEXT  J 

5725 

Tsat  =  <  T( 1 >+T<2) )*.5 

6730 

Tava=r<5> 
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o735  'jrac"37.'3853-*-.  :04388*7avg 

o/ziQ  Iarop"£nf  <  7)*1  ,  £+b/<  1  0  .  *Grad' 

8745  Tavac* T ( 5 ) -:aroo* . 5 

6750  IF  ABS<  Tav.g-7avgc>>z01  THEM 

6755  Tavg*< Tavg+Tavgc>* .5 

6760  GOTO  5705 

6765  END  IF 

§770- 

5775!  Compute  properties  of  'jater 

8730  Kw-FNKw(Tavg) 

5785  tluiua-FNMuw(Tavg) 

6790  Cpw  =  FNCpuj<  Tavg) 

6795  Prtd-FNPrtJ<Tavg) 

6800  Rhoui-FNRhoui(Tavg) 

5305! 

5310!  Compute   oroperties    of    Freon-114 

6815  Lmtd=Tdrop/L0G(<T(5)-Tsat)/<T(5)-Tdrop-T3at) ) 

6320  IF  Jj=0    THEN 

6325  Tu  =  Tsat+<rr*Lintd 

6330  ThetaD  =  Tij-Tsat 

6335  Jj  =  1 

6340  END  IF 

6845  Tf=(Tu+Tsat)*.5 

6350  Rho  =  frNRho(  Tf) 

6355  Mu=PNMu(Tf) 

6860  K-FNK(Tf) 

6865  Cp=FNCp<Tf) 

6370  3eta=^NBeta(Tf ) 

6875  Htg^NHfaUsat) 

5380  Ni-Ma/Rho 

6835  Alpna=</<Rho*Op> 

6390  Pr=Ni/Alpna 

6895! 

5900!  Analysis 

8905!  COMPUTE  MD0T 

6910  A=PI*<Do  2-Qi  '2)/4 

6915  P»PI*Do 

6320  rtdot=3.9657E-3+Fnis*<3.S1955E-3-Fms*<8.32006E-6-Fms*<  1  .23633E-7-Fms*^  .  31  397 

E-'0>  > ) 

8325  Q  =  1doc*Caiii*TaroD 

5920!  COMPUTE  NATURAL-CONNECTIVE  HEAT-TRANSFER  COEFFICIENT 

6925!  FOR  UNENHAMCED  END (3) 

6340  Hbar='90 

6345  Fe«<Hbar-P/<Kcu-A>> ".5*Lu 

6950  Tanh  =  frNTanh<Fe) 

5955  The ta= The tao* Tann/Fe 

6360  Xx-<9.31*Beta-Thetab*Do '3*Tanh/<Fe*Ni*Alpha>  >  '  . 166657 

6965  Yy  =  ( 1 -(  .559/Pr )    (3/V6> )    (3/27' 

6970  Hbarc  =  </Do*<  .6-*-.387*Xx/Yy )  '2 

6975  IF   ABS((Hbar-Hbarc)/Hbar)> .001    THEM 

6380  Hbar=<Hbar^nbarc>».5 

6985  GOTO  59*5 

5390  END  IF 

6995  ! 

7000!  COMPUTE  HEAT  LOSS  RATE  THROUGH  UNENHAMCED  ENDS 

7005  01  =<Hbar*p*t<cu.*A)  . 5*ThetaD*Tann 

7010  Qc=Q-2*Q1 

7015  As«*I»D2»L 

7020!  COMPUTE  ACTUAL  HEAT  FLUX 

7025  QdD=Qc/As 
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7020 
7035 
7040 

7045 

7050 

7055 

7060 

7065 

7070 

7075! 

7030 

7085 

7090 

7095 

7100 

7105 

7110 

7115 

7120 

7125 

7130 

7135 

7T40 

7145 

7150 

7155 

7150 

7165 

71  70 

7175 

7130 

7185 

7190 

7191 

7'95 

7200 

7205 

7210 

7215 

7220 

-»  -i  ~»  c 

I   _- J 

?220 
7235 
7240 
7245 
7250 
7255 
7260 
7255 
7270 
7275 
7230 
7?85 

7290 
7295 
7300 
7305 
7310 
7315 
7  320 


(1/3.  )*<Muua/FMMuu(Tui.)  ) 


IF    icfix-0    jR    Icfix-2    "HEN    Csf-!/Cf 

Thetab="Csf /Cp*Hf  g*<  Qdp/  <  Ma*Hf  q  >  ~<  .  0 ' 

Ho-Qap/  ThetaD 

0mega=Ho/Cf 

Uo=Q/<PI*Oo*L*Lmtd) 

Vw  =  1dot/<Rhouj*PI*Di  '2/4) 

Reut  =  ^houi»^ij'»D  i  /Muua 

Tui»Tw+0*Ruj/<PI*Oo*L) 

Gama=<u/Di-*Reu    .3*Pru 

PRINTER    IS    1 

Yij=(  1  .  /Uo-RuJ»0tnega 

Xu=Omeqa*Do/ <  Gama*0  i > 

Sx-Sx+Xu 

Sy  =  Sy  +  Yu 

Sxy  =  Sxy  +  Yu*Xuj 

Sx2-Sx2+Xui*Xu 

Sy2=Sy2+Yur*Yui 

NEXT    I 

ASSIGN    File    TO    * 

M=(Sx»Sy-Nru.n*Sxy)/(SxJ»Sx-Nru.TT*Sx2> 

C=<Sy-Sx-H)/Nrun 

IF  Icfix-0  OR  Icfix-3  OR  Icfix-4  THEM 

Cic-f/M 

Cfc-1/C 

END  IF 

IF  Icf Lx*1  THEM 

Cic-I/M 

Cfc-Cf 

END  IF 

IF  Icfix-2  THEM 

Cic«Ci 

Cfc-1/C 

END  IF 

IF  Icfix-3  THEN  7280 

IF  ABS( (Ci-Cic)/Cic)> 

Ci=<C.+Cic>-.5 

CF-<Cf+CFc)».5 

PRINTER  IS  1 

PRINT  USING  "2X 

PRINTER  IS  701 

GOTO  s6S5 

END  IF 

PRINT 

PRINTER  IS  701 

PRINT  USING  "23 

USING 

USING 


./Rr> 

(9.81 -Rho>  > 


5)  M/Rr)*P: 


001  OR  ABS((CF-Cfc)/Cfc)>  .001  THEN 


Csf  = 


:f 


MZ.3DE.2X 


,M.MZ.3DE":Csf 


2-2*1»Sxy-2-C-Sy+H"2*Sx; 


USING 
USING 


•10X. 

1  ox 


) 


PRINT 

PRINT 

PRINT 

5ur\2  = 

PRINT 

PRINT 

SUBEMD 

DEF  FMMuu< 

A=247.3/<  '  +  1 

Mu=2.4E-5*!0 

RETURN  *u 

-MEND 

DEF  FNCdi<j<T) 

Cpu-a. 211 20853- 

RETURN  Cpu*1 QUO 

FNEND 


'3X. ""ASSUMED       • 
•3X.  •"•CflLCULA  TED 


,MZ.3DE 
MZ.3DE 


'Mz. 


;DE":Cfa.C:a 
>DE";Csf ,Ci 


"Coet  f  t c len t 


■2*M»C*Sx+Nrun*C'*2 

>s    =    ""  .7  .  3DE"  :Sam2 

■    "••  ,D.DDD,':Rr 


•J  >J 

'A 


:5> 


'»/"->     "*  c  o  i  c  nr . 


-<  «4  .  -f< 


iEo+i 


J23E- 


T)  ) 
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7230  DEF  -"iRhou<  D 

7335  Ro-999.52946+T*<  .  0  !  263-"**  5  .  -4325  I  2E-  3-r*  I  .234147E-5)  ) 

7340  RETURN  Ro 

7345  FNEMD 

7250  DEF  FNPrum 

7355  Pru-FNCouj*  TWFNMuu»<T)/FMKu(T) 

7260  RETIJRM  Pru 

7365  FMEND 

7370  DEF  FNKui(T) 

7375  X  =  <T+272. !5)/273.  !5 

7230  Kw— .92247+X»<2.8395-X»<1 .300  7-X*<  . 52577- .  07344-X)  )  ) 

7335  RETIJRM  Kw 

7390  FMEND 

'0265  SUB  Select 
10370  C0,M  ,'Ido/'  lap 
10375  BEEP 
10330  PRINTER  iz    ■ 

10235  PRINT  USING  "4X .""Select  option: 

10290    PRINT    USING    "5X.""0    Taking    data   or    re-processing    previous    data' 

10335  PRINT  USING  ,,SX.',,,1  Plotting  data  on  Log-Log  

10400  PRINT  USING  "5X.',,,2  Plottma  data  on  Linear 

10405  PRINT  USING  "SX.""3  flaxe  cross-plot  coefft  file 

10410  PRINT  USING  "6X.""4  Re-circuiate  ojater 

10415  PRINT  USING  ,,5X.,,,,5  Purge 

10420  PRINT  USING  "6X."MS  T-Drop  correction 

10425  PRINT  USING  "SX,""7  Dnnt  Uo  rile 

10430  INPUT  IdP 
10435  IF  Idp-0  THEM  CALL  (lain 
10440  IF  Idp-1  THEN  CALL  Plot 
10445  IF  IdD  =  2  THEM  CALL  PI  in 
10450  IF  Iop  =  3  THEN  CALL  Coef 
10455  IF  Ido-4  THEM  CALL  flam 
10460  IF  IdP=5  THEN  CALL  3urg 
10465  IF  IdP=6  THEM  CALL  Tdcn 
10470  IF  IdP  =  7  THEN  CALi_  Uoprt 
10475  SUBEMD 


79 


APPENDIX  B 
FLOWMETER  CALIBRATION 

A  Fischer  Porter  Flowmeter  was  used  in  the  experimental 
apparatus  to  indicate  the  water  flowrate.  Prior  to 
conducting  the  boiling  tube  data  runs,  calibration  of  the 
flowmeter  was  performed  over  a  temperature  range  from  19  °C 
to  38  °C.  The  goals  of  the  calibration  were: 

1)  to  develop   a  mathematical  expression   relating  flowm- 
eter percent  reading  to  mass  flowrate  and 

2)  to   determine  a  viscosity   correction   factor  due   to 
varying  temperature  range  and  flowrates. 

A  weigh  tank,   platform- type  scale,   and  stopwatch  were  used 
to  record  water  collection  data. 

Data  were  processed  using  the  program  "FMCAL",  listed 
later  in  this  appendix.  This  program  took  the  flowmeter 
percent  reading,  water  weight  collected,  and  elapsed  water 
collection  time  as  the  inputs.  The  outputs  were  an  experi- 
mental mass  flowrate  and  a  mass  flowrate  difference.  This 
difference  was  given  by  a  mass  flowrate  computed  from  a 
correlation  minus  the  experimental  mass  flow  rate  for  the 
same  flowmeter  percent  reading.  The  correlation  was  based  on 
a  fourth-order  least-squares  fit  to  the  calibration  data. 
After  all  the  data  runs  were  completed,  the  data  points  were 
combined  into  one  file  and  a  single  fourth-order  polynomial 
was  generated  to  describe  the  relationship  between  the 
flowmeter  percent  reading  and  the  mass  flowrate.  Of  the 
original  data  points,  98.61  percent  were  within  a  ±5.5 
percent  range.  Reprocessing  the  data  points  led  to  78 
percent  within  a  ±  1.02  percent  range.  The  equation  (B. 1), 
resulting  from  the  calibration  data  analysis,  is  a  fourth- 
order  polynomial  from  a  curve  fit  of  the  reprocessed  data 
without  a  viscosity  correction  factor: 
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in  =  5.9637e-3  *    3  .  1690e-  3  (FMS)  (3.1) 


i 


-8.8020e-6(FMS)~  +  1 . 2369e- 7 ( FMS ) J 

-4 . 3190e-10(FMS)4 

(FMS  = flowmeter  percent  reading) 


Experimental  discharge  rates  are  less  than  the  rates 
calculated  from  theoretical  flow  equations  primarily  because 
of  the  effect  of  internal  molecular  friction  of  the  fluid 
and  the  viscosity.  The  theoretical  flow  equation  is: 


Q  =A..C 

»»  U 


w 


(3.2) 


where    Q=volumetric    flowrate,      A    =cross- sectional    area    of    the 

,v 

narrowest  part  of  annulus ,   V-=volume  of  float,   p c ~    density 

of  float,    X    =cross-sectional  area  of   the  largest   part  of 
i 

float,    and   c,.=density   of    fluid.    The   Coefficient   of 

w 

Discharge,  C-,  is  a  function  of  viscosity.   A  constant  value 
for  C^  in   equation  (3.2)   would  be  desirable   as  this  would 

GL 

show  a  negation  in  the  variations  due  to  the  viscosity.  The 
effect  of  the  float  shape  on  C  ,  with  varying  Re3.Ti.0lds  number 
is  shown  in  Figure  3.1  (a).  The  Reynolds  number  for  the 
flowmeter  uses  an  effective  diameter  of  the  tube  inner  diam- 
eter minus  the  float  outer  diameter.  This  plumb  bob  style 
float  presents  a  large  surface  area  to  the  fluid  stream  and 
the  viscous  effects  increase  as  the  flow  rate  increases. 
[Ref.  31:  pp.  9804-9808] 

The  data  collected  from  the  calibration  runs  showed 
little  or  no  change  in  the  Discharge  Coefficient  over  the  20 
to  90  percent  flowmeter  setting  range  and  the  15  K  change 
in  fluid  temperature.  The  square  edge  plumb  bob  float 
reaches   a  constant   Coefficient   of   Discharge  at   Reynolds 
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'igure  B.  1.  Effect  of  Float  Shaoe  on  Coefficient  of 

Discharge:  (a)  Snuare  Eaee  ?iumb  Bob  Float 
[trom  Ref.  30),  and  [b)    Taper  Edge  Plumb 
r loat. 
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number  >  10000.  This  is  due  to  the  effective  shape  of  the 
float  on  the  streamline  pattern  as  shown  in  Figures  B. 2  (a) 
and  (b).  The  experimental  calibration  range  of  1300  <  Re  < 
27000  displayed  a  flattening  out  of  the  Coefficient  of 
Discharge.  A  different  effective  shape  is  seen  by  the 
streamline  pattern  because  of  the  minimal  clearances  between 
the  tube  bead  guides  and  the  float.  Also,  the  float  used  in 
this  particular  flowmeter  had  tapered  edges  versus  square 
edges.  This  different  effective  shape  is  shown  Figures  B. 2 
(c)  and  (d).  Figure  B.  1  (b)  shows  the  effect  of  the  float 
shape  on  the  Coefficient  of  Discharge  for  the  tapered  edge 
float. 

Additionally,  the  accuracy  of  the  flowmeter  contributes 
to  the  flattening  effect  on  the  Coefficient  of  Discharge 
shown  by  the  experimental  data.  The  accuracy  effect  is  felt 
in  two  major  components,  reproducibility  and  scale  factor. 
Reproducibility  accuracy  is  minimally  determined  by  the 
change  in  the  fluid  flowrate,  corresponding  float  movement, 
and  observation  of  float  position.  The  calibration  data  had 
good  reproducibility  over  the  various  flowrate  and  tempera- 
ture ranges.  The  scale  factor  accuracy  is  inherent  to  the 
scale  markings  on  the  flow  tube.  Fischer  Porter  conducted 
numerous  tests  of  flowmeter  calibration  readings.  They  used 
a  maximum  scale  fraction  error  of  0.  9  mm  on  a  250  mm  tube. 
Based  on  their  results,  Fischer  Porter  reported  [Ref.  31: 
pp.  9814-9816]  an  accuracy  of  plus  or  minus  one  percent  at 
high  flow  rates  was  a  reasonable  expectation.  Also,  at  low 
flow  rates  ( <  20  percent)  a  resonable  expectation  would  be 
an  accuracy  of  eight  to  ten  perecent. 

The  experimental  data  followed  a  trend  of  negligible 
viscosity  effect.  The  flow  tube  design,  float  shape,  and 
scale  factor  accuracy  were  the  main  contributors  to  this 
negation  of  a  viscosity  correction  factor.  As  such,  equation 
(B. 1)  was  used  without  a  viscosity  correction  factor  to 
calculate  the  mass  flowrate  for  the  boiling  tube  runs. 
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Figure  B. 2.  Effect  of  Float  Shape  on  Streamline  Pattern: 

(a)  Square  cage  Fiumb  Bob  Float  (  trom  Ref.  30), 
and  ( b )    Taper^Edge  Plumb  Bob  Float. 
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APPENDIX  Q 
MODIFIED  WILSON  PLOT 

As  referenced  in  Chapters  III  and  IV,  this  appendix 
provides  the  data  reduction  procedures  used  to  compute 
inside  and  outside  heat- transfer  coefficients.  The  tubes 
used  for  data  collection  were  the  smooth,  modified  High 
Flux,  and  Modified  Turbo-3  tubes.  The  modified  Wilson  plot 
anaylsis  could  not  be  directly  run  on  the  enhanced  tubes 
since  reliable  boiling  correlations  for  these  surfaces  do 
not  exist.  Also,  only  two  unknowns  can  be  allowed  between 
correlations  for  the  inside  and  outside  heat- transfer  coef- 
ficients. The  steps  used  in  this  procedeure  are  outlined 
below  and  are  contained  in  DRP4,  •  which  is  cited  in  Appendix 
A: 

Basic  Equations 

The  overall   heat- transfer  coefficient   in  terms   of  the 
overall  thermal  resistance  is  given  by: 

:  J-%-     *     R    +—  (CI) 


'J 

o 


where 


(C.2) 


Modified  Wilson  Plot 


1.   Assume  C  r    and  Tf  for  the  boiling  side  and  compute: 


Z(  P-.-O 


0.  5 


n  = 


S,0! 


.-2 
uA  i 

.  2. 


(C.  3) 


Csfi/r 


(C.4) 


Note:  with  r  equal  to  0.333,  equation  (C.3)  represents  the 
Rohsenow  correlation  [Ref.  23:  p.  969] . 

2.   Assume  C  j_  and  compute: 

0  .1^ 
k      «  C .3  _  0 .23  /  u   1  ( r 


h .  =  c.r 
i    i 


(C.6) 


3.   S 


ubstitute  ho   and  h±      from  steps  1   and  2   above  into 
quation  (C.l)  and  rearrange  to  yield: 


eqi 


-  a 


a.  r 


+  r 


(C.  7) 


4.   Now,  let: 


Y  = 


w 


(C.8) 


and 


QA 


;<  = 


(C.9) 


Construct  the   least- squares  line   for  Y   versus  X   in  the 
form  of: 


Y  =  mX  *    C 


(C.  10) 
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5.   Compute  a  new  set  of  values  for  C .  and  C  -  as  follows 

1         Sx 


(C.ll) 


and 


sf 


(C.  12) 


6.  Repeat  steps  1  through  5  until  convergence  for  Ci  and 
Csf  between  two  successive  iterations  is  less  than  0.  1 
percent. 


APPENDIX  D 
DATA  RUNS 

The  following  table  outlines  the   data  runs  used  in  this 
investigation: 
Where: 

1.  WH  -  Smooth  tube 

2.  HF-  High  Flux  (i.e. ,  porous- coated  Korodense)  tube 

3.  HFM  -  Modified  High  Flux  (i.e. ,  porous  coating  machined 
off)  tube 

4.  TB  -  Turbo-B  tube 

5.  TBM   -  Modified   Turbo-B   (i.e. ,  external   enhancement 
machined  off)  tube 

6.  T  -  inlet  water  temperature 

7.  V  -  water  velocity 

8.  const  -  constant 

9.  dec  -  decreasing 
10.   inc  -  increasing 
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TABLE  4 
DATA  RUN  DESCRIPTION 


ile  Name 

Oil  % 

#  points 

Description 

WH01 

0 

10 

Light  Off  Effect 
( cold  water) 
Light  Off  Effect 

WH02 

0 

10 

(warm  water) 
Light  Off  Effect 

WHO  3 

0 

9 

(steam) 

const  V   dec 

WH05 

0 

10 

WH06 

0 

8 

T    const  V   dec 

WHO  7 

1 

10 

T    const  V   dec 

WH08 

1 

10 

T    const  V   dec 

WHO  9 

2 

10 

T    const  V   dec 

WHIO 

2 

10 

T    const  V   dec 

WHIOA 

2 

10 

Repeatability  WH10 

WH11 

6 

10 

const  V   dec 

WH12 

6 

8 

T    const  V   dec 

WH13 

6 

8 

Light  Off  Effect 
(warm  water) 

const  V   dec 

HF12 

0 

10 

HF13 

0 

12 

T    const  V   dec 
Light  Off  Effect 
( steam) 
Light  Off  Effect 

HF14 

0 

12 

( cold  water) 
Light  Off  Effect 

HF15 

0 

10 

(warm  water) 
Repeatability  of  1 

HF16 

0 

12 

1F12 

HF22 

0 

7 

T    inc  V   const 

HF23 

0 

6 

dec  V   const 

HF24 

1 

12 

const  V   dec 

HF25 

1 

10 

const  V   dec 

HF26 

1 

6 

dec  V   const 

HF27 

1 

7 

T    inc  V   const 

HF28 

2 

12 

const  V   dec 

HF29 

2 

12 

T    const  V   dec 

HF30 

2 

6 

T    dec  V   const 
(Repeatability  of 
T    inc  V   const 

HF2 

HF31 

2 

8 

HF32 

6 

10 

const  V   dec 

HF33 

6 

10 

const  V   dec 

HF34 

6 

7 

dec  V   const 

HF35 

6 

7 

T    inc  V   const 

HFMOl 

0 

8 

Sieder-Tate  Coeff 

Run 

HFM02 

0 

8 

Sieder-Tate  Coeff 

Run 

HFM0  3 

0 

10 

Sieder-Tate  Coeff 

Run 

HFM04 

0 

8 

Sieder-Tate  Coeff 

Run 

HFM05 

0 

8 

Sieder-Tate  Coeff 

Run 

HFM06 

0 

8 

Sieder-Tate  Coeff 

Run 

HFM0  7 

0 

10 

Sieder-Tate  Coeff 

Run 

HFM08 

0 

8 

Sieder-Tate  Coeff 

Run 

TBOl 

0 

10 

Light  Off  Effect 
( cold  water) 
Light  Off  Effect 

TB02 

0 

10 

(warm  water) 
Light  Off  Effect 

TB03 

0 

10 

( steam) 
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TB04       0       10        T    const  V   dec 

(^Repeatability  of  TB03) 
dec  V  const 
const  V  dec 
inc  V  const 
const  V  dec 
dec  V  const 
const  V  dec 
inc  V  const 
const  V  dec 
const  V  dec 
(Repeatability  of  TB12) 
dec  V  const 
const  V  dec 
inc  V  const 
const  V  dec 
dec  V  const 
inc  V  const 
const  V  dec 
jht  Off  Effect 
ream) 

TB21       6       10        Light  Off  Effect 


TB05 

0 

6 

t 

TB06 

0 

10 

T 

TB07 

0 

7 

T 

TB08 

1 

10 

T 

TB09 

1 

6 

T 

TB10 

1 

10 

T 

TB11 

1 

6 

T 

TB12 

2 

10 

"  T 

TB13 

2 

10 

T 

TB14 

2 

7 

TB15 

2 

10 

T 

TB16 

2 

6 

T 

TB17 

6 

10 

T 

TB18 

6 

7 

T 

TB19 

6 

7 

T 

TB20 

6 

8 

T 

ft 

TBM04 

0 

8 

TBM05 

0 

8 

TBM06 

0 

6 

TBM07 

0 

6 

TBM08 

0 

8 

TBM09 

0 

8 

TBM10 

0 

8 

( cold  waterl 


TB22       6       10        Light  Off  Effect 

(warm  water) 
Sieder-Tate  Coeff  Run 
Sieder-Tate  Coeff  Run 
Sieder-Tate  Coeff  Run 
Sieder-Tate  Coeff  Run 
Sieder-Tate  Coeff  Run 
Sieder-Tate  Coeff  Run 
Sieder-Tate  Coeff  Run 
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UNCERTAINITY  ANALYSIS 


The  uncertainity  for  mass  flow  race,  Reynolds  number, 
heat  flux,  LMTD ,  wall  resistance,  overall  heat- transfer 
coefficient,  inside  heat- transfer  coefficient,  and  outside 
heat- transfer  coefficient  were  analyzed  for  selected  run3  of 
the  smooth,  High  Flux  and  Turbo-3  tubes.  The  analysis  was 
based  on  the  Kline-McClintock  [Ref.  32:  p.  3]  method  of 
uncertainity  analysis.  For  example,  the  following  equations 
were  used  for  the  uncertainity  of  the  heat  flux: 


q  = 


q 


TIU_L. 


(S.1) 


and 


Q  s  Ac  CT.  - 

^     a   in 


OU" 


(2.2) 


In  accordance  with  Kline  and  McClintock,   the  uncertainities 
are  given  by: 


oq 


5L  Y 


0  .  5 


(E.  3) 


and 


50 


Q  |_\ 


*e 


./. 


1(7  -  r   ■ 

x   in    out 


57 


OU"! 


out 


0.5 


(E-4) 


The  uncertainity  in  the  conduction  losses  for  the  unenhanced 
ends  of  the  boiling  tube  were  considered  negligible  in 
comparison  to  the  boiling  surface  uncertainity  and  as  such, 
were  disregarded.  Table  5  lists  the  uncertainities  for  the 
previously  mentioned  values. 
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TA3LZ  5 
UNCERTAINTY  ANALYSIS  PERCENT S 

file     WHO  5     WHO  5     HF13     HF13      T304      T304 
q      31140    26690    49540    33780     63730     53670 
2-24     1.39     2.24     1.39      1.37      1.07 


(w/m    ) 
(m/s) 

3%e 
Re 
i3. 

sun 
utrn 

■4 
3 


1.56  2.51  1.56  2.51  1.43  2.51 

1.79  2.65  1.73  2.61  1.61  2.51 

1.59  2.52  1.59  2.52  1.46  2.52 

4.96  3.59  3.10  2.45  2.29  1.56 

2.76  2.75  11.17  11.17  2.92  2.92 


■■d 
f!k  5.21         '4.39  3.43  3.52  2.71  3.02 


h. 

i 

5h 

i 


2.03  2.56  1.61  2.25  1.43  2.21 

5.34  5.03  4.90  4.39  3.35  3.96 


Note:  All  runs  were  performed  in  purs  refrigerant 
WH  -    Smooch  cube 

H7  =  Koroaer.se  cube  with  High  Flux  Coating 
T3  =  Turbo -3  tube 

The  large  uncertainity  for  the  wall  resistance  of  file 
HF13  (porous-  coated  Korodense  tube)  is  due  to  the  10 
percent  uncertainity  of  the  thermal  conductivity  coeffi- 
cient. For  this  investigation,  the  coefficient  used  was  an 
average  value  for  Zhe  values  given  in  the  open  literature 
for  copper -nickel  olates  or  tubes. 
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APPENDIX  F 
LIST  OF  NOMENCLATURE 

NOMENCLATURE 

A  -  Surface  Area 

cn  -  Specific  heat  at  constant  pressure 

r. 

Cc  -  Oil  concentration 

C^  -  Inside  Sieder-Tate- type  coefficient 

C-.P  -  Rohsenow  coefficient 

bi 

d,-  -  Inside  diameter  of  tube 

d0  -  Diameter  of  mouth  of  reentrant  cavity 

D  -  Diameter 

e  -  Depth  of  internal  ridge  or  fin 

E,,  -  Energy  transfer  per  unit  volume 

f  -  Friction  factor 

g  -  Acceleration  due  to  gravity 

g^  -  Gravitational  constant 

h  -  Convective  heat- transfer  coefficient 

hfT-  Specific  enthalpy  of  vapor-liquid  mixture 

k  -  Conductive  heat- transfer  coefficient 

km  -  Conductive  heat  -  transfer  coefficient  of  metal 

L  -  Tube  length 

L0  -  Entrance  effect  length 

m  -  Mass  flow  rate 

p  -  Distance  between  ridge  or  fin  peaks 

Pr  -  Prandtl  number 

q  -  Heat  flux 

Q  -  Total  heat  transfer 

r  -  Tube  radius 

Re  -    Pveynolds    number 

Rw  -    Thermal    wall    resistance 

St  -    Stanton   number 
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uc 

V 

A 
a 
a 


P  " 
pr" 


Overall  heat- transfer  coefficient 

Specific  volume  of  vapor 

Velocity 

Change  (i.e.,  Tin  -  Tcut  ) 

Surface  tension 

Surface  tension  of  oil 

Surface  tension  of  refrigerant 

Viscosity 

Viscosity  of  oil 

Viscosity  of  refrigerant 

Density 

Density  of  oil 

Density  of  refrigerant 


SUBSCRIPTS 
f   -  fluid 
i   -  inside 
in  -  inlet 
m   -  mixture 

0  -  outside 

01  -  oil-liquid 
out-  outlet 

rl  -  refrigerant- liquid 

s   -  smooth  tube 

sat-  saturation 

w   -  water 

wo  -  outside  wall 


94 


LIST  OF  REFERENCES 


1.  Bergles,   A.   E.   and  Webb,   R.    L. ,    A  Guide  To  The 
Literature  on  Convective   Heat  Transfer  Augmentation, 
Advances  in   Enhanced  Heat   Transfer.   Vol.    43,   pp. 
81-89,  198b. 

2.  Wanniarachchi ,  A.  S.  ,  Marto,  P.   J.  ,  and  Reilly,  J. T.  , 

The  Effect  of  Oil  Contamination  on  the  Nucleate 
Pool-Boiling  Performance  of  R-114  from  a  Porous-Coated 
Surface,  ASHRAE  Transactions.  pp.  2-73,  1986 
(accepted  for  publication). 

3.  ASHRAE  Guide  and  Data  Book.  Chapter  18,   pp.   281-291, 
TW5T. 

4.  Dougherty,  R.  L.  and  Sauer,  Jr.  ,  H.  J.  ,  "Nucleate  Pool 
Boiling   of   Refrigerant-Oil    Mixtures   From   Tubes, 
ASHRAE  Transactions.  Vol.  80,  pp.  175-178,  1975. 

5.  Chongrungreong,  S.   and  Sauer  Jr.  ,   H.   J.  ,   "Nucleate 
Boiling     Performance     of     Refrigerants      and 
Refrigerant-Oil   Mixtures,"  ASME   Transactions.    Vol. 
102,  pp.  701-705,  November,  TWW. 

6.  Fujii,  M. ,  Nishiyama,  E. ,  and  Yamanaka,  G. ,   "Nucleate 
Pool  Boiling   Heat  Transfer  From  Micro-Porous  Heating 
Surfaces,   Advances   in  Enhanced   Heat  Transfer.    pp. 
45-50,  1§79: 


7.  Webb,  R.  L. ,  "The  Evolution  of  Enhanced  Surface 
Geometries  for  Nucleate  Boiling,"  Heat  Transfer 
Engineering.  Vol.  2,  nos  3-4,  pp.  46-69 ,  Jan- June , 
1981. 

8.  Reilly,  J.  T. ,  The  Influence  of  Oil  Contamination  on 
the  Nucleate  PooT-Boiling  Behavior  of  R-114  from  a 
Structured  Surface"  M~~  Si  Thesis  ,  Naval  Postgraduate 
School ,  Monterey,  California,  pp.  30-35,65,70,  March, 
1985 . 

9.  Czikk,  A.  M.  and  O'Neill,  P.  S.  ,  "Correlation  of 
Nucleate  Boiling  From  Porous  Metal  Films,"  Advances  in 
Enhanced  Heat  Transfer,  pp.  53-59,  1979. 

10.  Carnavos,  T.  C. ,  "An  Experimental  Study:  Pool  Boiling 
R-ll  With  Augmented  Tubes,  Advances  in  Enhanced  Heat 
Transfer.  pp7  103-108,  19$1. 

11.  Czikk,  A.  M. ,  Gottzmann,  C.  F. .  Ragi ,  E.   G. ,  Withers, 
J.   G. ,  and  Habdas ,   E.   P. ,    Performance  of  Advanced 
Heat   Transfer   Tubes   in   Refrigerant-Flooded   Liquid 
Coolers,"  ASHRAE  Transactions.  Vol.   76,   part  1,   pp. 
96-107,  1970": 


95 


12.  Arshad,   J.    and  Thome,    J.   R.  ,  'Enhanced  Boiling 
Surfaces:   Heat   Transfer  Mechanism  Mixture  Boiling, 
ASME-JSME  Thermal  Energy  Engineering  Joint  Conf. .  Vol. 
1,  pp.  191-197,  1983.    

13.  Arai ,  N.  ,  Nakajima,  T. ,  Fukushima,  T.  Fujie,  K.  Arai , 
A.  ,  and  Nakayama,  Y.  ,  "Heat  Transfer  Tubes  Enhancing 
Boiling  and  Condensation  in  Heat  Exchangers  of  a 
Refrigerating  Machine,"  ASHRAE  Transactions.  Vol.  83, 
part  2,  pp. 58-62,  197?. 

14.  Ayub,  Z.  H.  and  Bergles,  A.  E. ,  "Pool  Boiling  from 
Gewa  Surfaces  in  Water  and  R-113,"  Augmentation  of 
Heat  Transfer  in  Energy  Systems.  Vol.  52 ,  pp.  57-64 , 
1985 . 

15.  Nakayama,  W.  ,  Daikoku,  T.  ,  Kuwahara,  H.  ,  and  Nakajima, 
T.  ,  "Dynamic  Model  of  Enhanced  Boiling  Heat  Transfer 
on  Porous  Surfaces,"  Advances  in  Enhanced  Heat 
Transfer,  pp.  31-43,  1979. 

16.  Jensen,  M.  K.  and  Jackman,  D.  L.  ,  "Prediction  of 
Nucleate  Pool  Boiling  Heat  Transfer  Coefficients  of 
Refrigerant-Oil  Mixtures,"  ASME  Transactions.  Vol  106, 
pp.  184-190,  1984. 

17.  Marto,  P.  J.  ,  Reilly,  D.  J.  ,  and  Fenner,  J.  H.  ;  "An 
Experimental  Comparison  of  Enhanced  Heat  Transfer 
Condenser  Tubing, '  Advances  in  Enhanced  Heat  Transfer. 
pp. 1-8,  1979. 

18.  GA  Technoligies  Report,  GA-A17833,  Fluid  Mechanics  and 
Heat  Transfer  Spiral  Fluted  Tubing,  by  J.  S~i  Yampolsky 
eT.     al.  ,  1984. 

19.  DSR  Report,  70790-69,  Investigation  of  Heat  Transfer 
Augmentation  Through  ~Use  of  Internally  Finned  and 
Roughened  Tubes,  by  A.  E~.  Bergles  et~!  al.  ,  pp.  18-23 , 
197$, 

20.  Karasbun,  M.  ,  An  Experimental  Apparatus  to  Study 
Nucleate  Pool  Boiling  of  R-114  and  Oil  Mixtures.  MT 
S.  Thesis ,  Naval  Postgraduate  School ,  Monterey , 
California,  pp.  24-32,  54-56,  December,  1984. 

21.  Chaddock,  J.  B. ,  "Influence  of  Oil  Refrigerant 
Evaporator  Performance,'  ASHRAE  Transactions.  Vol.  82, 
part  1,  pp.  474-486,  1§76: 

22.  AERE  Report,  R7318,  The  Development  of  Heat  Transfer 
Tubes,  by  I.  H.  Newson  and  T.  D.     Hudgson,  19/3. 

23.  Rohsenow,  W.  M. ,  "A  Method  of  CorrelatingM  Heat 
Transfer  Data  for  Surface  Boiling  of  Liquids,"  ASME 
Transactions.  Vol.  74,  p. 969,  1952. 


96 


24.  Incropera,  F.  P.  and  DeWitt,  D.  P. ,  Fundamentals  of 
Heat  Transfer.  John  Wiley  and  Sons,  p.  ^06 ,  1981. 

25.  Withers,  J.  G. ,  "Tube-side  Heat  Transfer  and  Pressure 
Drop  for  Tubes  Having  Internal  Helical  Ridging  with 
Turbulent/Transitional   Flow   of   Single-phase   Fluid: 

Part    l--Single-helix    Ridging,"    Heat Transfer 

Engineering .  Vol  2,  No.  2,  1980. 

26.  Withers,  J.   G.  ,  "Tube-side  Heat  Transfer  and  Pressure 
Drop  for   Tubes  Having   Internal  Helical   Ridging  with 
Turbulent/Transitional   Flow   of   Single-phase   Fluid: 
Part    2--Multiple-helix   Ridging,     Heat    Transfer 
Engineering .  Vol  2,  No.  2,  198u 

27.  Nobukatsu,  A.  and  et.  al.  ,  "  Heat  Transfer  Tubes 
Enhancing  Boiling  and  Condensation  in  Heat  Exchangers 
of  a  Refrigerating  Machine,"  ASHRAE  Transactions.  Vol. 
83,  part  2?  pp.  58-70,  197?. 


28. 


Green,  G.  H.  ,  "Influence  of  Oil  on  Boiling  Heat 
Transfer  and  Pressure  Drop  in  Refrigerants  12  and  22," 
ASHRAE  Journal,  pp.  57-61,1982. 


29.  Stephan,  K. ,  Influence  of  Oil  on  Heat  Transfer  of 
Boiling  Freon  12  (Refrigerant  12)  and  Freon  22 
(Refrigerant  22),  Proceedings  of  Xlth  International 
Congress  of  Refrigeration.  Vol.  i,  pp!  369-3/9 ,  1963. 

30.  Stephan,  K. ,  and  Mitrovic,  J. ,  "Heat  Transfer  in 
Natural  Convection  Boiling  of  Refrigerant-Oil 
Mixtures,"  Models  and  Correlations  Review,  pp.  73-85, 
1982.       

31.  Fischer  and  Porter  Company,  Theory  of  the  Flowrator. 
Catalog  Section  98-A,  pp.  9801-9816,  194/. 

32.  Kline,   S.   J.   and  McClintock,   F.   A. ,   "Describing 
Uncertainities     in     Single-Sample     Experiments, 
Mechanical  Engineering,  p.  3,  Jan.  1953. 


97 


INITIAL  DISTRIBUTION  LIST 


No.   Copies 


1.  Library  Code  0142  2 
Naval  Postgraduate  School 

Monterey,  California  93943 

2.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  Virginia  22314-6145 

3.  Department  Chairman,  Code  6  9  1 
Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943 

4.  Professor  Paul  J.  Marto,  Code  69  Mx  3 
Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943 

5.  Dr.  A.  S.  Wanniarachchi ,  Code  69Wa  2 
Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943 

6.  Dr.  B.  H.  Hwang,  Code  2722V  1 
David  W.  Taylor  Naval  Ship  Research 

and  Development  Center 
Annapolis,  Maryland  21402 

7.  Mr.  R.  Helmick,  Code  2745  1 
David  W.  Taylor  Naval  Ship  Research 

and  Development  Center 
Annapolis,  Maryland  21402 

8.  Mr.  N.  R.  Clevinger  1 
Product  and  Process  Development  Mgr. 

Wolverine  Tube 
P.  0.  Box  2202 
Decatur,  Alabama  35602 

9.  Mr.  J.  McManus  1 
Box  139 

Logan  Lake,  British  Columbia 
Canada  V0K  1W0 

10.  Lt.  S.  McManus  2 
6840  N.  Columbus  Blvd 

Tucson,  Arizona  85718 

11.  Mr.  Elias  Ragi  1 
Union  Carbide  Corporation 

P.  0.  Box  44 

Tonawanda,  New  York  14151 

12.  Dr.  Masafumi  Katsuta  1 
School  of  Science  and  Engineering 

3-4-1,  Ohkubo,  Shinjukuku, 
Tokyo,  160,  Japan 


98 

1  7  8  8  8 


ALP0::  CHO( 


2l95«i9 


^Thesis 
&2537 


HcManus       i-boilitvg 
nucleate  p       t  and 

°£  *"  wtures  from  water- 
oil  mixtures       £aceS. 
heated  enhanced 


2195*t9 


Thesis 

M2537 

c.l 


McMarms 

Nucleate  pool-boiling 
of  R-114  refrigerant  and 
oil  mixtures  from  water- 
heated  enhanced  surfaces, 


